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Abstract 

The totally asymmetric simple exclusion process (TASEP) on Z 
with the Bernoulli-/} measure as initial conditions, < p < 1, is sta- 
tionary. It is known that along the characteristic line, the current 
fluctuates as of order t 1 / 3 . The limiting distribution has also been ob- 
tained explicitly. In this paper we determine the limiting multi-point 
distribution of the current fluctuations moving away from the charac- 
teristics by the order t 2 ^ 3 . The main tool is the analysis of a related 
directed last percolation model. We also discuss the process limit in 
tandem queues in equilibrium. 

1 Introduction and result 

Continuous time TASEP. The totally asymmetric simple exclusion process 
(TASEP) is the simplest non-reversible interacting stochastic particle sys- 
tem. In TASEP, particles are on the lattice of integers, Z, with at most one 
particle at each site (exclusion principle). The dynamics is defined as fol- 
lows. Particles jump to the neighboring right site with rate 1 provided that 



*Department of Mathematics, University of Michigan, Ann Arbor, MI, 48109, USA 
e-mail: baik@umich.edu 

^University of Bonn, Endenicher Allee 60, 53115 Bonn, Germany 
e-mail: ferrari@uni-bonn.de 

*Institut Fourier, 100 Rue des maths, 38402 Saint Martin d'Heres, France 
e-mail: Sandrine .Peche@ujf-grenoble . f r 



1 



the site is empty. Jumps are independent of each other and take place after 
an exponential waiting time with mean 1, which is counted from the time 
instant when the right neighbor site is empty. 

It is known that the only translation invariant stationary measures are 
Bernoulli product measures with a given density p G [0,1] (see [21] )■ In 
the sequel we fix a p G (0, 1) to avoid the trivial cases p = (no particles) 
and p = 1 (all sites occupied). One quantity of interest is the fluctuations 
of the currents of particles during a large time t. Consider the so-called 
characteristic line given by x = (1 — 2p)t. Then, the current fluctuation 
seen from the characteristic are (9(t 1//3 ) (in agreement with the scaling for 
the two-point function established in [5]), and the limit distribution function 
has been determined (conjectured in [2S] and proved in [T7], see also [T]). 
The limiting distribution was first discovered in the context of a directed last 
passage percolation model [I]. 

On the other hand, if one looks at the current along a line different from 
the characteristic line, say x = ct with c / 1 — 2p, then for large time t one 
just sees the Gaussian fluctuations from the initial condition on a t 1 ! 2 scale, 
since the dynamics generates fluctuations only (^(t 1 / 3 ), thus irrelevant [14]. 

A non-trivial interplay between the dynamically generated fluctuations 
and the one in the initial condition occurs in a region of order t 2 / 3 around the 
characteristic line. One of the main results of this paper is the determination 
of multi-point limits of the current fluctuations around the characteristics 
x = (1 - 2p)t + 0(t 2/3 ) (see Theorem O] and Theorem [T7]). 

There has been a great deal of work pertaining to the limiting fluctuation 
of TASEP over the last ten years since the well-known work of Johansson [T9] . 
See for example, the review paper [16] for the limit processes which arise from 
deterministic initial conditions. A recent paper [5J, building on the earlier 
work [10], is concerned on the situation where random and deterministic 
initial condition are both present, but not stationary. 

Directed percolation model. It is well-known that the currents of TASEP 
(of arbitrary initial condition) can be expressible in terms of the last passage 
time of an associated directed last passage percolation model (see for exam- 
ple, P2J,28J). We first discuss the asymptotic result of the following directed 
percolation model associated to stationary TASEP J2S]: see the paragraphs 
preceding Theorem 11.61 below for the exact relation to stationary TASEP. Let 
w i,ji hj — 0, i,j G Z, be independent random variables with the following 
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distributions 

^0,0 = 0, 

w it0 ~ Exp (1/(1- p)), i>l, , , 

w od ~Exp(l/p), j>l, 1 • ] 

Wij ~ Exp(l), i,j > 1. 

Here the notation X ~ Exp(r) means that X is a random variable expo- 
nentially distributed with expectation r. An up-right path ir from (0, 0) 
to (x, y) G N 2 is a sequence of points G Z 2 , 1 = 0, . . . , x + y), starting 
from the origin, 7r = (0,0), ending at (x,y), ir x + y = (x,y), and satisfying 
n £+ i - n £ e {(1,0), (0, 1)}. Denote by L(ir) = Yj(i,j)e% w iJ- Then > the last 
passage time is defined by 

G(x,y) = max L(n). (1.2) 

We are interested in the limit distribution of the properly rescaled last pas- 
sage time of G([xN], [yN]) in the N — > oo limit. 

When there are no 'borders' i.e. w^j = when % = or j = 0, the limiting 
distribution of G(x,y) was first obtained in [19]. The case of 'single-border' 
when u>j )0 = but Wqj = Exp(l/p) was studied in [2]. The above 'double- 
border' case was considered in [T7(|28] . (The Poisson and geometric variations 
were studied earlier in jl].) One can also consider a more general model when 
w ifi ~ Exp(l/(l/2 + a)) for % > 1 and w 0J ~ Exp(l/(l/2 + 6)) for j > 1, 
where a,b > —1/2. Such model was considered in [2B] and a conjecture 
on the limiting distribution left in that paper was recently confirmed in pQ. 
For the Poisson and geometric variations, this generalization was considered 
in [3] earlier when x = y. 

An interesting characteristic behavior in such bordered models is the tran- 
sition phenomenon. One can imagine that the last passage time comes from 
the competition between the contributions from the 'bulk' i,j > 0, and the 
edges i = or j = 0. See for example, Section 6 of [2] for an illustration 
of such heuristic ideas. For the model (11.11) . a crucial role is played by the 
critical direction (which corresponds to the characteristic line of TASEP), 

x (1 — p) 2 ^ 

It is easy to check that along a direction other than the critical direction, the 
fluctuations of G([xN], [yN]) is given by Gaussian distribution on the N l l 2 
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scale, see Appendix |D] This situation corresponds to the Gaussian fluctua- 
tions along a non-characteristic line in a stationary TASEP [J3] mentioned 
earlier. 

Along the critical direction and also in a N 2 ^ 3 neighborhood of the critical 
direction, the fluctuations of G([xN], [yN]) are on the N 1 ^ 3 scale. Precisely, 
the following, among other things, is proven in [T7]. Let us set 



x ■■= p(i - p)- 



;i-4) 



Then set the parameter N 
Consider the scaling 



[(1 — ^X)T\ , where T is considered to be large. 



x(r) 

y( r ) 

(r, s) = T - t 



l_ p )2 T + r j^l r 2/3 



1-2% 

p 2 T _ r iA^ T 2/3 



1-2* 
2(1 - 2p) X 1/3 
l-2x 



;i.5) 



rp2/3 _|_ 



T l/3 



A 



1/3- 



The parameter r measures the displacement of the focus with respect to the 
critical line (on a T 2 ! 3 scale). In particular, for r = one looks exactly along 
the critical direction. £(t, s = 0) is the macroscopic value of the last passage 
time for large T, while the parameter s in £{t, s) measures the amount of the 
fluctuations (on a T 1 / 3 scale) of the last passage time. Then for any fixed 
TGl, 

lim F{G{x(r),y(r)) < £(r,s)) = F T (s) (1.6) 

T->oo 

for an explicit distribution function F T , which satisfies F T (s) = F- r (s). An 
analogous result was first obtained in [3] for a Poissonized version of the 
problem in which F T was obtained in terms of a solution to the Painleve II 
equation (see (3.22) of [4J: F T {x) = H(s + r 2 ; r, — r)). Based on this result, 
Prahofer and Spohn conjectured in [2H] that (11.61) holds. This conjecture was 
proven in [T7| whose analysis was somewhat different from [3] . This results in 
a different formula of F T , expressed in terms of the Airy function (see (1.20) 
in [IT]). It can be checked that these two formulas do agree. 

One of the main objects of study in this paper is the limit of multi-point 
distribution G(xj,yj), j = 1, • ■ ■ ,k. The limit of the process for the no- 
border case was first obtained in [21] following the earlier work of [21] on the 
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Poissonized version of the model. The limit of the process for the double- 
border case was considered in [IT] for the model when w ,o ~ Exp(l/(a + b)) 
andwj )0 ~ Exp(l/(l/2+a)) fori > 1, w 0J ~ Exp(l/(l/2 + 6)) for j > 1 when 
(a, b) e (-1/2, 1/2) and a + b > 0. Note that when a = -b = 1/2 - p, the 
random variable wq^ becomes singular in this model: the restriction u>o,o = 
is significant in connection to stationary TASEP. 

The geometric counterpart of the model (II. ip was studied earlier by Ima- 
mura and Sasamoto [18]. Denoting by Geom(g) a random variable with the 
probability mass function (1 — q)q k , k = 0,1,2, ■■■ , the authors of [18] consid- 
ered the model (11.11) where u7 ,o = 0, u^o ~ Geom(7 + a), woj ~ Geom(7_a) 
and Wij ~ Geom(a 2 ) for i, j > 1. Since the exponential model (II .ip can be 
obtained as a limiting case of a — > 1, the analysis of [TB] can in principle be 
used to yield the corresponding result for the model (11.11) . Nevertheless, this 
paper differs from [18J on the following aspects, (a) The authors of [18] ob- 
tained explicit limiting distribution functions for the case that the situation 
a + b > (and a + b < 0) in (II. ip . However, they left the case corresponding 
to a + b = as the limit of the a + b > case and did not compute the limiting 
distribution explicitly (see remarks after Theorem 5.1 in [IB])- This critical 
case is the most interesting (and the most difficult) case for our situation 
and we give an explicit formula in Theorem 11.21 below, (b) The justification 
of the limits of Fredholm determinant and other quantities appearing in the 
analysis requires proper conjugations of corresponding operators in order to 
make sense of the Fredholm determinant and trace class limit. These issues 
were not discussed in [18] (the main issue was to determine the possible limit 
regimes and not specifically the a + b = case), (c) In addition to the multi- 
point distribution on the line x + y — constant considered in [18], we also 
obtain limit of the process result for points (xk,Vk) not necessarily on the 
same line (see Theorem II .5p . (d) The limit of the process for points at more 
general positions than on a line mentioned in (c) is used to prove the limit 
of the process (in the sense of finite distribution) of stationary TASEP (see 
Theorem 11.61 and Theorem 11.71 below) . 

We first state the result extending (I1.6P to the joint distributions at points 
on the line 

C N :={(x,y)>0\x + y = N} (1.7) 
at and near the critical direction. We first need some definitions. 
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Definition 1.1. Fix m G N. For real numbers T\ < r 2 < . . . < r m and 

^1) • • • > "5m; 5Ct 



oo /-oo 

2^ (»+!/) 



7^ = si + e 3 r i / rf x / dy Ai(x + y + if )e~ 
dx Ai(ar + y + Tf)e~ TiX , 

POO POO 

= e -K dX dye- x{Tl - n) e- Tiy Ai{x + r f 2 + A) Ai(y + r\ + A) 

JO J sy 

2_ f*S\ — 2 POO 

+ 1^2] i-l^= = / dye'^D- dy Ai(y + x + if )e 7iV . 

A/47r(Tj - n) J-oo JO 

(1-8) 

/or z, j = 1,2, ... ,m, where Ai denotes the Airy function. 

The first result is the following. 

Theorem 1.2. Fza; m G N. For real numbers T\ < r 2 < . . . < r m and 
Sx, • • • , s m , with the scaling given in (ll.5p . 

Km p( f){G(x(r fe ),y(r fc )) <£(r fc ,s fe )} j 

^ fc=1 ' (1.9) 

= y2l^[9 m (r, s) det (l - P s £ Ai P s ) J • 

t^<9Sfc V V /L 2 ({l,...,m}xR)/ 



fc=l 

"2/ 



i/ere L ({1, . . . , m} x M) is equipped with the standard measure v®dx where 
v is the counting measure on {1, . . . .m}. P s denotes the projection operator 
P s {k,x) = t[ x>Sk ], and Kai is the so-called extended Airy kernel JW$ with 
shifted entries defined by the kernel 

K A i((i,x),(j,y)) := [K M ]ij(x,y) 

■ POO 

J d\Ai(x + \ + T?)Ai(y + \ + Tf)e- x ^- T *\ ifn^Tj, (L1Q) 
- f d\ Ai(x + A + r?) Ai{y + A + Tf)e~ x ^- Ti \ if n > Tj . 

J — oo 

The function g m {r, s) is defined by 
g m (r, s) = K+(pP s $,P s y) 

m m POO POO /I 1 1 \ 

= U + J2J2 dx dy^ j (y)p 3ii (y,x)^ i (x), 

l=\ , = 1 Js i Js 3 



where 



p:= {t-P s K M P s y\ Pj,i(y,x) :=p((j,y),(i,x)), (1.12) 

and $((i,x)) := $i(x), ^((j, y)) = ^j(y)- Finally the functions 1Z, $ ; and 
\I/ are defined in Definition \l.l\ 

Observe that dist((a;(r/ c ), y(r k )), C N ) < 2. 

Remark 1.3. When m — 1, fll.9p agrees with the limiting function in (1.20) 
of [T7], as one may expect. 

Remark 1.4. The fact that 1 — P s Kx\P s is invertible follows from the fact 
that P s Ka[P s is trace class (see [21]) and that det(l — P S K^P S ) > for all 
given s6i. See Lemma [B. II in Appendix iBl 

The shift in the integrand by rf is due to the fact that the last passage 
time is (macroscopically) a linear function along the line £jy, in contrast 
with the non-border case (i.e., Wi^ = u> j = 0) where the last passage time 
has a non-zero curvature. Therefore, when \rk\ ^> 1, the contribution from 
det (l-PK Ai P) will very close to one and the main contribution comes from 

9m(T, S). 

The second theorem is a generalization of Theorem 11.21 to the case when 
the points (xk,yk) are not necessarily on the same line Cn, x + y = N. 
We show that the fluctuation is unchanged even if some of the points are 
away from the line to the order smaller than 0(T). This is due to the slow 
de-correlation phenomena obtained in |15j : along the critical direction the 
fluctuations decorrelate to order C(T 1//3 ) over a time scale 0(T) (instead of 
C(T 2/3 )). More precisely, if we compare the last passage time G at two points 
(x, y) and (x', y') with (x' — x,y' — y) = r ■ ((1 — p) 2 , p 2 ), their fluctuation will 
be r + 0{r 1 ^) (see Lemma [5.31 below for details). Consider a v e (0, 1) and 
any fixed number 9 and consider the scaling (a generalization of (ll.5jl ) given 
by 



x{t,6) 



y(r, 



(1 - p) 2 (T + 6T V ) + r ^—T 2 ^ 

1 - 2 X 

p 2 (T + 9T V ) - t 2x4/3 T 2 / 3 
P K ' 1 - 2x 



;i.i3) 



2fl-2oW 1 / 3 , T 1 / 3 
£(r, 8,s) = T + 6T U — r Z[ \ Zp J X T 2 / 3 + s —tr- 

See Figure 11.11 for an illustration. 
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Figure 1.1: Assume that the black dots are 0(T U ) for some v < 1 away from 
the line Cn- Then, the fluctuations of the passage time at the locations of 
the black dots are, on the T 1 / 3 scale, the same as the one of their projection 
along the critical direction to the line C^, the white dots. 



Theorem 1.5. Fix m e N and v e (0, 1). For real numbers T\ < r 2 < . . . < 
r m , 9\, . . . , 9 m , and si, . . . , s m , with the scaling given in f ll . 13|) . 

Km P t f]{G(x(r k , 9 k ), y{r k , 9 k )) < £{r k , 9 k , s k )}\ 

\k=i J (114) 

m d ( \ 
= Y / ^(9m(T,s)det(l-P s K Ai P s ) ). 

9s k V V / L2({l,...,m}xR)y 

This generalization of Theorem 11.21 is proven in Section [5j 

Stationary TASEP and directed percolation. We now discuss the result in 
terms of stationary TASEP. The mapping between TASEP and last passage 
percolation model is as follows. We assign label to the particle sitting at 
the smallest positive integer site initially. For the rest we use the right-to-left 
ordering so that • ■ ■ < x 2 (0) < xi(0) < < x (0) < x_i(0) < • ■ ■ . Then 
Xfc(t) > Xfc + i(i) for all t > 0, since the TASEP preserves the ordering of the 
particles. 

For i,j 6 Z such that i — j > Xj(0), let j) be the waiting time of 
particle with label j to jump from site i —j — 1 to site i —j (the waiting time is 
counted from the instant where particle can jump, i.e., particle is at i — j — 1 
and site i — j is empty). The j) are iid Exp(l) random variables. Let 
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L(x, y) be the last passage time to (x, y) E N 2 along a directed path in the 
domain D := E < x,j < y,i — j > Xj(0)}, starting any point in 

the domain. This is a curve-to-point optimization problem. An example of 
the non-trivial part of the boundary of this domain is illustrated in Figure [L2l 
below as the curve at t — 0. Then a straightforward generalization of the 
step-initial case of [TH] shows that 

P(fT =1 {L(z fc>3/fc ) < t k }) = n^ii^vM >x k - y k }). (1.15) 

When the initial condition is random, the domain of directed percolation 
is also random. Note that (see Figure [L2|) . the part of the domain D in the 
first quadrant is the rectangle {1 < i < x, 1 < j < y}, but the parts in the 
second quadrant % < 0,j > 1 and the forth quadrant i > 1, j < are of 
random shape. Observe that D does not intersect with the third quadrant 
i,j<0. 

Define — (C- + 1) to be the right-most empty site in {. . . , —2, —1} in 
the initial particles' configuration and to be the position of the left-most 
particle in {0,1,...}. First let us focus on the case (_ = 0, ( + = (as 
in Figure II. 2p . Then, since the initial condition is stationary Bernoulli, 
an interpretation of Burke's theorem [12], see also [13], shows that the 
{L(0, j)\l < j < y} is distributed as {X u Xi + X 2 , . . . , X l + ■ ■ ■ + X y } 
where Xj's are iid Exp(l/p) distributed. Similarly, by considering holes 
instead of particles one finds that {L(z,0)|l < i < x} is distributed as 
{Yi, Y\ + Y 2 , . . . , Y\ + ■ ■ ■ + Y X } where l^'s are iid Exp(l/(1 — p)) random vari- 
ables. Hence L(x,y) has the same distribution as G(x,y) defined from ( II. ip . 
This argument was sketched in Section 2 of [28]. Consequently, we have for 
x k ,Vk > 1, h > 0. 

n^T=i{^y k (h)>x k -y k })=n^T=i{G(x k ,y k )<t k }), (1.16) 

for stationary TASEP when initially the position is empty and the position 1 
is occupied. 

For generic £_ > and ( + > 0, L(x,y) is distributed as G(x,y) where 
now we assume in ( II. ip that w^q = ■ ■ ■ = w^ +t0 = iy ,i = • • • = w ^_ = 
where ( + ~ Geom(l — p) and ~ Geom(p). But as shown in Proposition 
2.2 in [17], this change does not affect the asymptotics. 

Geometrically, the TASEP and directed percolation can be thought of as 
two different cuts of the three-dimensional object 

{x,y,G(x,y)\x,y>l}: (1.17) 
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Characteristics 



Figure 1.2: The dots are (random) particle configurations at time t — and 
some later time t. The position of particle with label n = y is its projection on 
the J-axis. Interpolating between particles as in this example, one gets a line 
configuration, which is interpreted as height function ht(j) above position j. 



(a) the directed percolation problem we analyze is the cut at {x + y = N}, 

(b) the particles' configuration of the TASEP at time t is the cut at {G = t}. 

We have shown in Theorem 1 1 . 51 that the limit process does not depend on the 
cut chosen for the analysis as long as we avoid the cut along the characteristic, 
{x/y = (1 — p) 2 / p 2 }. Thus, to get the fluctuations around the characteristic 
line it is enough to project on {x + y = (1 — 2x)t} (see Figure 11.21 for an 
illustration) , for which the limit theorem was proven in Theorem 11.21 

The following two theorems are proven in Section [5j Define the functions 

n{r)= [p 2 T-2rp X 1/3 T 2 / 3 \, 

q{r) = L(l - P fT + 2r X V3 T 2/3 _ (1 _ p)sT ^/ x ^\ . (L18) 

Theorem 1.6 (Particles' position representation). Fix m G N. For real 
numbers T\ < r 2 < . . . < r m and s 1; . . . , s m , 




\ m d 

n(r k )(T) > q{r k )} J = (dUr, s) det (l - P s K Ai P s 

1.19) 



ds k 
k=i K 
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An equivalent but geometrically slightly different way of representing the 
TASEP is via a height function (see e.g. [T?ll2"g] ). Let us define the occupation 
variable, f]i(t) = 1 if there is a particle at site i at time t and rji(t) = 
otherwise. Then, define the height function 

f 2JV t + £Li(l-27*(t)) J for j>l, 
h t (j) = I 2N t , for j = 0, (1.20) 

{ W-YLj+ii 1 -**®)' forj<-l, 

where N t is the number of particles which jumped from site to site 1 
during the time-span [Q,t\. Then the link between the height functions and 
the locations of particles is 

nn™ =1 {h tk (x k - yk) >x k + y k }) = P(nr =1 {x, fc (t fe ) > x k - y k }). (1.21) 

Let us consider the scaling 

J(r) = L(l-2p)T + 2r X 1 /3 r 2/3j ) 

H{r) = - 2 X )T + 2r(l - 2p) X 1/3 T 2 / 3 - 2 SX 2/3 T 1/3 J . 

Theorem 1.7 (Height function representation). Fix m e N. For rea/ nam- 
6ers ri < r 2 < . . . < r m and sj_, . . . , s m . Then, 

tm \ to, q 

f]{h T (J(r k )) > H(r k )} = ^ — (<? m (<r, s) det (l - P a K M P a ty . 
k=l / k=l k 

(1.23) 

Remark 1.8. For simplicity, in Theorems 1 1 . 6 1 and 1 1 . 71 we stated the result 
only for fixed time. However, the statements can be extended to different 
times in a similar manner of the extension from Theorem 11.21 to Theorem 1 1.5 1 



Queues in tandem. There is a direct relation between queues in tan- 
dem and TASEP. Suppose that there are infinitely many servers, and we 
assume that the service time of a customer at each server is independent and 
distributed as Exp(l). Once a customer is served at the server i, then the 
customer joins at the (z+l)th queue, and so on. In other words, the departure 
process from the ith queue is the arrival process at the (i + l)th queue. Sup- 
pose that the system is in equilibrium with parameter p: the arrival process 
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Figure 1.3: Queues in tandem in equilibrium. The black dots represents the 
costumers and at every white dots one changes to the next counter. 



at each queue is independent Poisson process of rate p. Then the departure 
process at each queue is also independent Poisson process of rate p due to 
Burke's theorem (see for example, [121E3]; see also [25]). This also implies 
that at each time, the number of customers in each queue is distributed as 
Geom*(l - p). (Here X ~ Geom*(l - p) means that ¥(X = k) = p(l - p) k , 
k = 0,1,2, ■■■ .) Now consider a fixed time t = and arbitrary select one cus- 
tomer and assign label to that customer. For convenience, we call the queue 
in which that customer is in at time as the Oth queue. We assign labels to 
the other customers so that the labels decreases for the customers ahead in 
the queues (see Figure [L~3]) . Let Qj(t) denote the label of the queue in which 
the jth customer is in at time t. The mapping from the queueing model to 
the TASEP is obtained by setting Xj(£) = Qj(t) — j (see Figure fl~3l . The 
equilibrium condition implies that the initial condition for the corresponding 
TASEP is stationary. Hence if we define Ej(i) be the time the jth customer 
exits the queue i, then we find that 

F(nU{E yk (x k - i) < t k }) = nnUiQy k (tk) > x k }) 

= nn™ =1 {*y k (t k )>x k -y k }). 
Hence using (11.161) and Theorem 11.51 we immediately obtain the following 



;i.24) 
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result. 



Theorem 1.9. Fix m G N and v G (0, 1). For real numbers T\ < T2 < . . . < 
r m , 9i, . . . , 9 m , and s±, . . . , s m , with the scaling given in fll . 13|) . 



To prove Theorem 11.21 which is the basis for Theorems 11.61 and 11.71 we first 
consider a slightly different directed percolation model which is known to 
be determinantal. This model is then related to ( II. ip by a shift argument 
(Section [2]) followed by an analytic continuation (Section [3]), resulting in 
an explicit formula of the joint distribution of the last passage times of a 
more general version of (II. ip when w i0 ~ Exp(l/(l/2 + a)) for % > 1 and 



w 0J ~ Exp(l/(l/2 + b)) for j > 1, for a,b G (-1/2, 1/2) (see Theorem E£2J. 



The formula for the model (II .ip is obtained by setting a = — b = p — 1/2. 
These arguments are a generalization of the arguments given in [4] and [17] 
for one-point distribution. In Section H] we carry out the asymptotic analy- 
sis of the formula obtained in Theorem 13.21 The proofs of Theorems ll.2| 11.61 
and II .71 are given in Section |5j Some technical computations are given in Ap- 
pendices: various expressions of integrals involving Airy functions are given 
in Appendix |Aj the invertibility of an operator appearing in Theorem 11.21 is 
discussed in Appendix [Bj certain operator are shown to be trace-class in Ap- 
pendix and finally we explain the Gaussian fluctuation along non-critical 
directions in Appendix [Dj 

Remark 1.10. As mentioned above, Theorem 13.21 below contains a for- 
mula for the model with u>j j0 ~ Exp(l/(l/2 + a)) for i > 1 and Wqj ~ 
Exp(l/(l/2 + b)) for j > 1, for a,b G (—1/2,1/2), which is more general 
than (II. ip . This would correspond to other random initial data which possi- 
bly allows a shock. Essentially all the ingredients for the asymptotic analysis 
for this more general case are in this paper, and one can obtain the limit 
laws for such models. However, for the sake of simplicity, we do not pursue 
this direction here. 




(1.25) 



Outline 
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2 The shift argument 

We first consider a slightly different directed percolation model. Let 

w 0fi ~ Exp(l/(a + b)), 

uJ i)0 ~Exp(l/(l/2 + 6)), i>l, 

w ,;~Exp(l/(l/2 + a)), j > 1, 1 ' > 

w id ~ Exp(l), i,j > 1, 

where the parameters a and b satisfy 

a, b G (-1/2,1/2), a + b>0. (2.2) 

Denote by G^ b (x,y) to be the last passage time from (0,0) to (x,y) for 
this modified model. This model is well-studied and has nice mathematical 
structure. In particular, the correlation functions and joint distribution 
functions on are determinantafj], which is well-suited for an asymptotic 
analysis. Note that the original model given in (11 .ip corresponds to the case 
when a + b = and {y ,o — 0. We will show in Sections [2] and [3] how to obtain 
a joint distribution formula for the original model ( II .ip from this modified 
model ( 12. ip . Let G a> t,(x, y) be the last passage time for the model ( 12 .ip with 
wofi replaced by 0. We proceed as follows: 

(1) Shift argument: for a, b G (—1/2, 1/2) with a + b > 0, we relate the 
distribution of G a ^ with the one of G^ b . 

(2) Analytic continuation: we determine an expression for G a ^ which 
can be analytically continued in all a,b G (—1/2, 1/2). 

(3) Choice of parameter: finally we set a = p — 1/2 and b = 1/2 — p. 

1 It is a limit of geometric random variables studied in [201 HI] or of the Schur mea- 
sure 26 . Exponential random variables are studied directly in [11] on a more general case 
than the present one. 
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The step (1) is done in Section [2] and step (2) is presented in Section [31 

Proposition 2.1 (Shift argument). Let a,b G (—1/2,1/2) with a + b > 0. 
Let (xx,yi), . . . , (x m ,y m ) be the a set distinct of points in I? + and define 

/ m > 
P(«l, . . .,U m ) ■= P ( P|{G a)b (x fe ,2/ fc ) < Mfc} 

(2.3) 

P + ( Ul , . . . , u m ) := P ( (^{G^Or*, y*) <"*})■ 



vfc=l 



(1 m <9 A 
i+^E^J^K- ■■■«-)■ <") 

Proof of Proposition [2~7ll Let us only consider m = 2: the proof for general 
m > 2 is a straightforward generalization and we leave it to the reader. Set 
r = a + b. Then 

POO 

P + (u 1 ,u 2 ) = / d2/P(G+ 6 (xi,yi) < u u G+ b (x 2 ,y 2 ) < u 2 \w 0fi = y)f{w 0fi = y) 
Jo 

POO 

= / dy P(ui- y,u 2 - y)re~ ry . 
Jo 

(2.5) 

Consider the Laplace transform 



du x I du 2 P + ( Ul ,u 2 )e- tlUl - t2U2 
Jo 

POO POO poo 

r dy dux du 2 P( Ul -y,u 2 -y)e- ry - tlUl - t2U2 (2.6) 
Jo Jo Jo 

poo poo poo 

r dy dz x \ dz 2 P(z u z 2 )e- ry - tl{zi+y) - t2iz2+lj) . 



jv J -y J -y 

Since P(z\, z 2 ) = when either z\ < or z 2 < 0, we can restrict the integral 
from —y < Zi < oo to < Zi < oo. Then the above equals 



OO POO 



dy I dz x I dz 2 P( Zl ,z 2 )e-( r+tl+t2)y ~ tlZl - t2Z2 

poo (2.7) 
dz x I dz 2 P(z 1 ,z 2 )e~ tlZl - t2Z2 . 



r 



r + ti + t 2 Jo 
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Multiplying by t±h±h anc j integrating by parts leads to 

/ du x / du 2 P(u u u 2 )e- tlUl - t2U2 
Jo Jo 

rOO rOO 

/ du x / du 2 P + ( Ul ,u 2 )e- tlUl - t2U2 
Jo Jo 



r + tx + t 2 



dux J™ du 2 (l + t - 1 ^ r ^j P + (ux, u 2 )e- hu i- t2U * 

oo poo / I f d d \ 
dux / du 2 P + (ux, u 2 ) + - 7; h ^ — P + (ux, u 2 ) 1 e" 



o Jo 



r V dux du 2 



t\U\—t2U2 



(2-8) 

Since this holds for all t 1 ,t 2 > 0, by inverting the Laplace transform, we 
obtain 

P( Ul ,u 2 ) = P + ( Ul ,u 2 ) + ±(J- + ^)P + (ux,u 2 ). (2.9) 

□ 

As mentioned earlier, the probability P + (ui, . . . ,u m ) has an explicit de- 
terminantal expression. In fact, the joint distribution function at points 
on Cn for the directed percolation model with w^j ~ Exp(l/(a» + bj)), 
di + bj > 0, has a determinantal structure (it is a limit of geometric random 
variables studied in [20j[2T] or of the Schur measure [26]; see [7] for a direct 
approach to exponential random variables) Specifically, this follows, for ex- 
ample, from Theorem 3.14 and (1.3) of [21 J after substituting aj i— y 1 — ^ and 
bi i — y 1 — in (1.3) and taking the limit L — > oo. The kernel in this limit is 
explicitly derived in Theorem 3 of [UJ. For good survey papers on the topic, 
sec [22, 32J . By specializing to the case with weights as in ( 12. ip . we have the 
following result. Consider a set of distinct points (xx, yi), ■ ■ ■ , (x m , y m ) on the 
line L 2 t = G Z 2 | i + j = 2t}, which can be ordered and parametrized 

by tx < ••• <t m G [-*,*], 

x k = t + t k , y k = t-t k . (2.10) 

Set 



a+w) 



t-ti 



,2 ~ W 



2 The result can be extended to any set of points which can be connected by down-right 
paths, called space-like path, but we do not enter in the details here. See 7,8 for some 
examples in similar situations. 
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and[f| 

m = zl / tM£Z, gM i / dz (2 . 12) 

2m Jr a , 1/2 a-w 27ri/ r ^ i/2 z + 6 
Define the kernels 

K ij{. x >y) = Ki,j( x iy) - v i,j( x iv) ( 2 - 13 ) 

where 



2tti J iM (f)j[z) 
K id (x, y) = <£ dzl dw e zy ~ lr ! °' ( "' j 



(2.14) 



(2vri) 2 7 r _ £ ' w - z' 



-1/2 ^ 1 1/2 

Then, 

P + (n 1; . . . ,« m ) = det (1 - P u KP u ) LH{1 _ m}xW+) (2.15) 
where P u (k, x) = l[ x >u k ] and the operator K is defined by the kernel 

K((i,x),(j,y)) =K id (x,y) := K i>:j (x,y) + (a + b)f i (x)g j (y). (2.16) 

In (12.15p . we have used a slight abuse of notation as P U KP U is not a trace- 
class operator: one can indeed observe that for i > j, P Ui Vi t jP Uj {x,x) -/¥ 
as x — > oo. Nevertheless with a suitable multiplication operator M, the 
conjugate operator MP U KP U M~ X becomes trace-class so that the Fredholm 
determinant is well-defined, and the identity becomes valid analytically. More 
concretely, if we take M as in (I2.18P below, then it is shown in Appendix ICl 
below that MP U KP U M~ X is trace-class operator for a, b e (0,1/2). Thus, 
after conjugation, the Fredholm determinant is well-defined, and the identity 
becomes valid analytically, as proved in the following. 

Proposition 2.2. Let a, b G (0, 1/2). Fix constants a±, . . . , a m satisfying 

— a < a.\ < a 2 < ■ ■ ■ < a m < b. (2-17) 

Define the conjugated operator K con i by the kernel 

K%*(x,y) = M l (x)K hJ (x,y)M J (y)-\ M t {x) := e~ a > x . (2.18) 



3 For any set of points S, the notation ^ dzf(z) means that the integration path goes 
anticlockwise around the points in S but does not include any other poles. 
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Then, P u K con iP u is a trace-class operator on L 2 ({1, . . . ,m} x R) and 

P + ( Ul , ...,u m )= det (1 - PuK^P u ) L2({1 _ m}><R+) . (2.19) 

Proposition 12.21 and (I2.15P give the standard extension from one-point 
kernel to the multi-point (or extended) kernel. All the ingredients are in- 
cluded in Appendix B and Section 3 of [17], and also in [3~T] . The extended 
kernel can be found in the more recent work [llj. There the setting is more 
general allowing several lines/columns to have waiting times different from 
1. 

Observe K is independent of a and b, and the only dependence on a and b 
in K is through the rank-one term (a + b)fi(x)gj(y). Using this, we can find 
an expression of P + (ui, . . . , u m ) in which the condition (I2.17P can be relaxed 
to (I2.20p . This relaxation is important when we take the limit a + b — > in 
the next section. 

Proposition 2.3. Let a,b G (0, 1/2). Fix constants a%, . . . , a m satisfying 

1 1 , 

--<«!< a 2 <■■■ < a m < -. (2.20) 

Define the conjugated operator K c ° ai by the kernel 

K™?(x,y) = M i (x)K i>j (x,y)M j (y)~\ M^x) := e"^. (2.21) 

Then 

P + (u 1 , ...,u m ) 

= (1 - (a + b)((l - P u KP u y l P u f, P u g)) ■ det (l - PJ?^ P u ) ■ 

where the notation (, ) denotes the real inner product in L 2 ({1, . . . ,m} xl + ). 

Remark 2.4. One can check that the expression f)2.22p can be analyti- 
cally extended to a, b G (—1/2,1/2) with a + b > 0. Nevertheless since 
we will discuss the issue of extending the domain analyticity of a formula of 
P(ui, . . . ,u m ) to a, b G (—1/2, 1/2) (with no restriction of a + b > 0) in the 
next section, we do not discuss the detail here. 
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Proof of Proposition 12.31 From (12.161) . 



P u K co ^P u = Pjr° n3 P u + (a + 6)(P u / conj ) <g> {g^Pu) (2.23) 

where f con i and g con i are multiplication operators by the functions 
r^(i,x) = Mi{x)fi{x), g COQ Kj,y) = 9 S {V) M M~ X - The Unctions f coa i(i,x) 
and g con ^{j,y) are L 2 (M): see the end of the proof of Proposition IC.ll 

For the convenience of notations, set [K] = P u K con ^P u , [K] = P u K c ° n} P u , 
[f] = P u f coni and [g] = g co ^P u . Thus 

det (1 - [K]) = det (l - (a + b)(l - [K])"^/} <8> [g]) ■ det (l - [K]) , 
= (1 - (a + &)<(! - mr'ifl [9] )) • det (l - [K]) . 

(2.24) 

The above step holds assuming that 1 — [K] = 1 — P U K J P U is invertible. 
The invertibility can be verified as follows. Consider the modification of 
the directed percolation model where uiij = if i and/or j are equal to 
0, and Wij ~ Exp(l) for i, j > 1 (i.e., model without sources). Denote 
by P the new measure. This is the model with a = b = in (12. ip . and 
hence P(ui, . . . ,u m ) = det(l — P u K c ° ni P u ) for any given ux,...,u m > 0. 
It is easy to obtain a lower bound P{u\, . . . ,u m ) > (1 — e~ £ ) 2t > 0, with 
£ = min{«i, . . . , u m }/2t. Indeed, it is enough to take the configurations with 
Wij < minj-ux, . . . , u m }/2t for i, j > 1 such that i+j < 2t. This implies that 

det(l — P u K con} P u ) 7^ 0, and hence 1 — P u K c ° ni P u is invertible. 

Finally, since P U KP U is a bounded operator in L 2 (R), and 
P u f,P u g G L 2 (M), ((1 - \K])- l [fl [g] ) equals ((1 - P u KP u y l P u f,P u g) by 
conjugating back M. Now the condition (12. 17ft for the a^s can be relaxed 
to the condition ( 12.201) since K c ° ni is trace-class under this assumption: see 
Proposition IC.ll □ 



3 Analytic continuation 

We now find a formula extending (12. 4ft to the case when a,b G (—1/2, 1/2) 
without the restriction that a + b > 0. For this purpose, we show that both 
sides of ( 12. 4p are analytic in the parameter a, b G ( — 1/2,1/2). Analyticity 
of the left-hand-side of (12 .4p i.e., of P(ui, . . . ,u m ) is a straightforward gen- 
eralization of Proposition 5.1 in [T7]. This section is devoted to finding the 
analytic continuation of the right-hand-side of (12 ,4p . 
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By using (12.221) . the right-hand-side of (12. 4p becomes 



(a + 6) + ]T — J (-—-((i-P u KP u )-ip u f t p u9 )^j det (l-PjT^ P u ) . 

(3.1) 

The last determinant is independent of a, b. It is enough show that 

' ((l-P u KP u )- l P u f,P u g) (3.2) 



a + b 

is analytically continued in a, b G (—1/2,1/2). Note that by changing the 
contour, 

fi{x) = f^(x) + f^(x) (3.3) 

where 

= /f /2, M = ^/ (3.4) 

27T1 / r , a — w 

J 1 1/2 

Proposition 3.1 (Analytic continuation). For a,b G (0,1/2), 
1 



a + 6 

where 



({l-PuKP^PJ.P^) = -((l-PuKP^PvFtP^+Rajb. (3.5) 



P 1 / W 0i(q) e -u l( a- W ) 

«a,6 = 7T^ f dw j~T~W V HP 3 - 6 

27ri /r aj _ 6 ,_ 1/2 0i H (w - a) (w + 6) 



and F((z,x)) = -Fj(x) wift 

/*0O 

F l {x) = f^ /2 \x) + dyK hl (x,y)fi a \y) + l [l > 2] dyV iA (x,y)fi a \y). 

J U\ J ~oo 

(3-7) 

TTie term R a ,b is analytic in a, 6 E (—1/2, 1/2), and 

m /*OG 

((l-P u KP u )- l P u F,P u g) :=J2 dx ((1 - P^P^ 1 P u F)((i,x)) gi (x) 

i=l ^ 

(3.8) 

zs convergent and is analytic in a,b G (—1/2, 1/2). 

Hence the right-hand-side of (I3.5P is an analytic continuation of (13. 2p 
to a, b G (— 1,|). Combining (12 ,4p . (13. ip and (I3.5p . we finally obtain the 
following representation of P(ui, . . . , u n ) defined in (12. 3p . 
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Theorem 3.2. Recall the conditions and definitions from ( 12.1 Op 
through flZJl) a&oue. For a,b E (-1/2,1/2), 

P(«i, . . . ,u r 



5 

m 



a + 6 + 2 du~ k J [i Ra > b ~ - PuKPuy'PuF, P u g)) det (l - PjT° nj P, 

(3.9) 

for Mi, ... , u m G M+, where R a ,b and F are defined in Proposition ^. 11 

Proof of Proposition 13.11 Pari P Decomposition. First we decompose the 
contribution coming from the pole at 1/2 and a of /, namely 

m = ^Tb ~ ((1 -PuKPur'Puf^lPug) 

-{{t-PuKP^PJ^.P^). 

From Lemma E31 below, PJ^ = (1 + P U VP U )P U F^ + P U V(1 - P„)P (a) , 
where P^ is defined in (13TTS]) . Here, (1/(1 - P„)P( a ))(z, x) is well-defined 
pointwise and is in L 2 (R) as we can check as in the proof of Lemma 13.31 
Hence using the identity (recall that K — K — V in (j2.13j) ) 

(1 - P u KP u y\t + P U VP U ) = 1 + (1 - P u KP u y l P u KP u1 (3.11) 

the last term in ( 13 . 101) becomes 

(P u F^,P u g) + ((1 - P U KP U )-\P U KP U + P U V(1 - P u ))F^,P u g). (3.12) 

Observe that the function (KP U + V(t — P u ))F^ a \(i, x)) is precisely the last 
two terms in (13. 7p . Hence from (I3.10p . we obtain 

= -((1 - PjCP^P^P^) + — j- - (P u F^,P u g). (3.13) 

a + o 



Now a direct computation shows that 

POO 

(P u F^,P u g)= / dx f[ a \x) 9l {x) 

J ui 



2vri /r„ 1/2 ,_ 6 ^ + & JU1 
dw- 



(3.14) 



2?ri /r_ 1/2 ,_ 6 (w + b)(a-w) 

1 / , 0i(a) e -«i(«-«0 1 
4 dw ) ' - — r + 



2ni Jr-i,2,- b ,a <Pi(w) (w + b)(a - w) a + b' 
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Here in the third equality, we used the fact that the contour T_i/ 2 -b can 
be made to be on the left of the point w = a since a, b > 0. Hence ( 13. 13)) 
and fl3TT4j) imply fl33|) . 

Part II: Analyticity. We now show that the functions on the right-hand- 
side of (13.51) are analytic in a, b G (—§,§)• Clearly, i? a ,f> is analytic in 
a, b G (—1/2, 1/2) since both the poles w = a and w = —b lie inside the 
integration contour. We need to show that ((1 — P U KP U )~ 1 P U F, P u g) is ana- 
lytic. Note that K, K and V are independent of a, b. As f\ (x) = <pi(a)e~ ax 
is analytic in a G (— |, |), Fi(x) is analytic in a in the same domain. Also 
it is clear from the integral representation (12.121) that gi(y) is analytic in 
b G (— f> i)- Hence it it enough to show that ((1 — P U KP U )~ 1 P U F, P u g) is 
well-defined for a, b G (— |, |). Fix <5 £ (0, |). Let a, 6 G [— | + 5 , | — VI- 
Using the identity 

(1 - P u KP u y l P u F = (1 + P u ifP u (l - P u KP u y l )P u F, (3.15) 
and the estimates ( I3.23f) and ( I3.25f) in Lemma 13.41 below, we see that 

|((1 -PjTP u )-^P u F)(z,x)| < Ce-^ 2 -^', x > Ml (3.16) 

for some constant C > 0. On the other hand, from (13.241) . 

\gM\<Ce (1/2 ~ SM , yeR, (3.17) 

for some constant C > 0. Therefore, the inner product is convergent, and 
the Proposition is obtained. □ 

Lemma 3.3. Let a G (0, |). Define the function F^ in L 2 ({1, . . . , m} x M.) 
by 

Ft\x) = f[ a \x)5 iA . (3.18) 

Then VF^ defined by 

(VF^)((i,x))= [ dyV hl (x,y)fi a \y) (3.19) 



is well-defined for each x G IR and is in L 2 {M). Moreover, 

f(a) = F (a) +VF (a)_ ( 3 2 0) 
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Proof of Lemma D Recall from that fj a \x) = e~ ax (j)i{a). The esti- 
mate (I3.22p shows that the integral in (I3.19P is well-defined. (VF^)((i,x)) 
is well-defined pointwise. Note that Vi t i(x,y) = v(x — y) for some function 
v G L 2 (M)nL 1 (R) (see (109]) ). Hence the integral in (13191) equals {v*f[ a) ){x). 
Hence the L 2 (R) norm of the integral is bounded by H^Hl^R) II A IU 2 (R)- But 
as |v(a;)| < Ce-^-^M by f l3T22|) . we find that the integral in is in 



Now let r\\ G (0, a) and r/ 2 G (a, |) be fixed. In the integral formula of 
i(a;, 2/) in (I2.14p . we can change the contour iR to either iR + 7/1 or iR + 772. 
We will use the contour iR + r\\ when y > and the contour iR + 772 when 
y < 0. Then for i > 2, 

2m \J iR+m a-z(px{z) J m+m a-z(px{z) 
= e- ax ( j> l (a) = ft\x) 

(3.21) 

where the integral is evaluated using Cauchy's formula. Taking into account 
the case when i — 1, we obtain (13.201) . □ 

Lemma 3.4. For any 5 G (0, ~] ; there exists a constant C > such that 

\Vij(x,y)\ < Ce^ 2 - s ^-y\l [l>jh x,yeR, (3.22) 

and 

\Kij(x,y)\ < Ce- (l ' 2 -^ x+ y\ x, y > 0. (3.23) 
For any b G (— |, |) ; there is a constant C > swc/i i/iai 

|^(y)|<Ce-^ y G R. (3.24) 

Finally, suppose that a G (— |, |) zs given. For any 6 G (0, (| + a)(| — a)], 
there is a constant C > snc/i t/iat 

\Fi{x)\ < Ce-( 1/2 ~ S)X , x> Ui . (3.25) 

Here the constants are uniform if the parameters a, b, 5 are in compact sets. 

Proof of Lemma 13.41 When y — x > 0, by deforming the contour to 
iR - 1/2 + 5 in tfLWi . we obtain for i > j, 



\ Vi j(x,y)\=e-W-s)(y-*)±. 

< Ce -(l/2-S)(y-x) 



-1/2 + 5 + is) 



(3.26) 
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since the integrand is absolutely convergent (recall that —tj>l for i > j). 
We obtain the bound ( I3.22p similarly by using the contour iR + 1/2 — 5 when 
y — x < 0. 

The estimate (13.231) is easily obtained by taking the contours r_]/ 2 and 
Yi/2 as the circles of radii of 5 centered at —1/2 and 1/2, respectively, 

in dnsD. 

Recalling gj(y) (from (I2.12p ) we obtain the estimate (I3.24p by evaluating 
the residue at z = — b and making the remaining contour r_i/2 small enough. 

Finally, in order to estimate Fi(x), first note that \ f i \x)\ < Ce^^l 2 ^^ 
by using the integral representation (13 .4p with the contour given by 
the circle of radius 5 centered at 1/2. Also from (I3.22p and (I3.23p . 
\K i;1 (x,y)\ < C e -(i/2-5)(^)_ Using t h ese estimates and |A (a) (^)| < Ce~ ax 
in the definition (13. 7p . 

poo pUl 

F l (x) = ri 1/2 \x)+ / dyK hl (x,y)fi a \y) + l {l > 2] / dyV^{x,y)f[ a \y), 



(3.27) 

we obtain ( 13. 25 ft for x bounded below. In particular, to have the integrals 
over y bounded we need < 5 < | + a for the integral with and 
< 5 < \ — a for the integral with V^\. Both conditions are satisfied for 
0<5< (±-a)(± + a). □ 



4 Asymptotic analysis 

By setting a = — b = p — 1/2 in Theorem 13.21 we obtain 
P(ui, ...,u n ) 

(Ra,- a - ((1 - PuKPu^PuF, P u g)) det (l - PtX^Pu) ■ 

k=l " 

(4.1) 

We now begin asymptotic analysis of this formula. 

To obtain our main theorem (Theorem ll.2p we need to consider the fol- 
lowing scaling limit. Fix p e (0,1). Set x — p(l — P)i b — 1/2 — P, and 
a = p — 1/2. For a large parameter T, according to ( 11.51) and 1 12. lOft . we 




24 



consider 

_ 1 - 

' ~ _ 2 2 ' l-2x 



t = — T, U = -T + T t — T 2/3 , 



2fl - 2/jW 1 / 3 , T 1 / 3 ( ^ 4 ' 2 ' ) 
1 _ 2 X x 7 

where we order ri < r 2 < . . . < r m with t( G [— i, t] for all i. The convergence 
of the Fredholm determinants is ensured only after (yet another) proper 
conjugation. For this purpose, set 

A{i) = Z(i) exp(2r 3 /3 + r iSi ), Z(i) = ^(^e"^. (4.3) 

Proposition 4.1. Consider the scaling fiJJty . Then 



fT\ 173 A(j) ^ 

i im _ -Tr\ K hj( u i> u j) = [-^Ai]ij(sj,Sj), (4.4) 



uniformly for Sj, Sj in a bounded set. The operator is defined in U.10\) . 

Proof of Proposition 14.11 Recall that K iy j = K iy j — Vij, with V^j = for 
Ti < Tj. The same structure holds for Kai- Indeed, using the identity ( 1A.7I) 
of Lemma IA.ll we can rewrite ( 11.101) as follows, 



dX Ai(si + A + if) M( Sj + A + r?)e 
oxp \ + |(r| - rf ) + (r, Sj - r,*)) 



\(Tj—Ti) 



l{Ti > r 3 ) 

(4.5) 

The proof is divided into the convergence for tyy in Lemma 14.21 and of K^j 
in Lemma 14.41 below. □ 



Lemma 4.2. Consider the scaling h4-<fy an d i> 3- Then 

T\ 173 A(j) 

■ -\ ',.,!>/,. // ,) 



exp ( + + + ^ 



(4.6) 



47r(ri - 73) 



uniform for Si — s,- in a bounded set. 
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Proof of Lemma 14.21 Recall that Tj > Tj. We derive the asymptotics by 
saddle point analysis. Set 

g Q {z) := (r, - [(1 " 2p)z - vln(| - i 2 )] , 

ft (2) := -zfa - Si )v- 1/3 , (47) 

<?(z) :=<7oC*)+T- 1/3 <7i(4 
Then by plugging (0J2D into fl2Tl4j) . 

V id (ui,Uj) = — / dFexp (T 2/3 5-(F)) 

R (4 8) 

= -L f &p (T^g,{I) + TV* gi (z)) . 

There is a unique critical point for go in the interval (—1/2, 1/2), namely 

z c = a = p-1/2. (4.9) 
A straightforward computation gives 

^ = (n-r^(^ 75 L_ + (I75 L_), ,4.10) 

from which <7o(3" c ) = 2(7$ — Tj)Y~ 2 / 3 > 0. Also observe that 

e T 2 / 3 9 (a) = e T 2 /3g o(a)+T l/3 ffl(o) = ZQ_^ ^ n ^ 

For the saddle point analysis we use the contour 7 := {z c + it\t e R}. 
First let us show that the contribution coming from \t\ > 5 > is negligible 
in the T — y 00 limit. Denote by 7^ := {z G 7 | |Im(z)| > 5}. Then we have 



It'-- 



ft©" 3 ™) 

1/3 



- ^exp(Re[T 2 /^(i c + it)-T 2 /^(a)]) 



y\ 1 poo J 

) - dt- 



(4.12) 



X ; vr 7 5 (|1 - it/(l/2 - i" c )||l + it/(l/2 + 



T y/s i 



X7 ttA (l + t 2 /p 2 )^ 2 (l+t 2 /(l-p) 2 )"/ 2 
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with r\ := ti—tj = (ri — Tj) i^x ^ 2 ^ 3 — ^ wrien T is large enough. For p > 1/2, 
we have the bound 

1 1 

(l+t 2 /p 2 )"/2(l + t2/( 1 _ p )2^/2 ^ (\ ^' 13 ^ 

Hence using the linear lower bound 

1 + (t/pf > (1 + 5 2 /p 2 )(l + (t - 6)25/(5 2 + p 2 )) (4.14) 
for all t, 5, we have 

, *(TT7W- (1 + * 2/ ' ,2r, |^- (4 - 15) 

Therefore there exist a constant p = p(5) > and a constant C = C(p) > 
such that 

^C ^-ffP - (4,6) 

When p < 1/2, we obtain the same estimate by just replacing in some of the 
bounds p by 1 — p. 

Next we determine the contribution from a (^-neighborhood of the critical 
point. Noting the Taylor series 



go(z) = g (z c ) + (n - T^x-'^iz-ZcYil + 0{z-z c )), 
we find that 



9\{z) = gi{z c ) - (si - Sj)x 1/3 (z - z c ), 



(4.17) 



g(z) - g(z c ) - (r, - r,)x- 2/3 (?- ? c ) 2 + ( * " c) 



< sup 

teB(z c ,\z-z c \) 



«b w (t) 



Z - Zc 



3! 

/3 

3!(1-^ X ) V(l-P- l^-^cl) 3 ' 



4X 4/3 / 1 1 . , , , 

< ( r * - T j)—, — ( ~ 1 1~ + 77 i~ ~ in 3 ) l z " ^ 



(4.18) 

Thus by choosing 5 small enough, we find that for \z — z c \ < 5, 

dUBJ < ^^f\z-z c \ 2 , (4.19) 
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and also 

(jUBD < C\z-z,f. (4.20) 
for some constant C > 0. Using (I4.19p . (I4.20p and the general identity 

\e z - e w \ <\z-w\ max{|e z |, |e*"|}, (4.21) 

we deduce that 

/ ^ 2/3 ^) _ _L / 
'iWs 2n J-siT/x) 1 / 3 



i 

2^ 



"5(T/x) 1/3 

dF fe T2/3(9{a+iy(T/xrl/3) - 9(a)) - e -Y 2 (n-T 3 )-iY( Sl - Sj )\ 

-,i(T/x)V3 V / 
r-<5(27x) 1/3 

27rTV^ J_ s(T/x)l/3 

(4.22) 

Finally, since 

1 ^/X) 1/3 , 1 / (o._o.\2 

J- / ^v2c,.^.\_iv/..^,.'i ... -L / l* 5 — o ■ ' 



e ' J ('»- r i)- i n»i-«j)dy exp 



2VT J_ S( T/ x )i/3 x/47r(Ti - Tj) V 4 ( r i- T i) 



2tt J-SiT/x) 1 / 3 2l1 

< e- cTV3 

(4.23) 

for some constant c > 0, which is uniform for Sj — Sj in a bounded set, we 
obtain flU} from fT4TT6]) and (Q2"]l . □ 

Before discussing the asymptotic of K we obtain also an exponential 
bound of V it j(ui, Uj) when s», Sj — >■ +oo, which will be used later. 

Lemma 4.3. Consider the scaling \4-<fy an d i > j (i-e., Ti > Tj). Then, for 
any given k > 0, there exists a Tq > large enough such that 



T\ 1/3 A(j) 



< (J e -(n-r j )(s l +s j )/2 e - K \s l -s j \ (4.24) 



for all Si,Sj G R and T > T . The constant C is uniform in Si,Sj and in 
T>T . 
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Proof of Lemma 14.31 Set 



a + 



(*i-Sj)X 1/3 
2(T i -T j )TV3' 



a+ ^ 3 

a+ 2(r<-T J )> 



a 



1/3 



if I Si - Sj | < eT 1 / 3 , 
if Si - sj > eT 1 / 3 , 
if Si - Sj < eT 1 / 3 , 



(4.25) 



where e > is a fixed constant chosen small enough so that z c lies 
in a compact subset of (—1/2,1/2). Taking the integration path as 
7 := {z c + it 1 1 G R}, a computation as in f)4.12p yields (the function g(z) 
is defined in ( I4.7p ) 



Z{i) 



^ J dlexp (T 2 / 3 ^ - T 2 ' 3 g{a)) 



< J_ / (// . 

2tt 



(4.26) 



(|l-i*/(l/2-* c )||l + i*/(l/2 + z c )|)* 
where 77 = (7* — T j)^z^T 2 ^ 3 > 1 when T is large enough. Hence 

fOHn< fi ^/3 g ( Ze) - T2 /3 g(a) i r <// 

7T y 



D T 2 /3 9 ( 2c )-T 2 /3 ff (a) 



1 / 1 



The last integral satisfies 



dt 



7T V 2 



1 + 



o (l + tV(l/2-|^i)^ 

2 d 



eft 



(4.27) 



(1 + * 2 K 



ds 



e s ds 



(4.28) 



as r] = 0(T 2 / 3 ) — > 00. Therefore we find that there is a constant C > such 
that for T large enough, 



T 



\ i/3 



Z(j) 

xj m 



- 



<C-e W (T 2 / 3 (g(z c )-g(a))) 



(4.29) 



Now from (14.181) and (I4.19p (with z = z c and z c = a), if we have taken e 
small enough, then 



g(zc)-g(a)< 2 ^ 2/3 (z c -a) - xl/3rV3 



^ - a 



(4.30) 
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Now we plugging the value of z c in (I4.25P for three difference cases. When 
\si-sA <eT l '\ 



g(z c ) - g{a) < 



[Si - Sj, 



[Ti ~ Tj 



)T 2 / 3 ' 



(4.31) 



When Si - sj > eT 1/3 , then 



g(z c ) - g(a) < 



e(A( Si - sj) - ZeT 1 / 3 ) 



< 



e{Si - s j/ 



.( r ._ r .) T l/3 - 8 ( r ._ ri ) T l/3- 



(4.32) 



When Si - Sj < -eT 1/3 , then 



g(z c ) - g(a) < 

Therefore, we obtain 
T\ 1/3 Z(j) 



E(4\8i-Sj\ -3eT^ 3 ) 



< 



£ | S>i S j | 



\(Ti - - 8(r t - Tj )TV 3 



(4.33) 



X J Z (i) Vi ^ Ui,Uj/ 



^•exp(-|^), \ Si - Sj \<eT^, 



< 



C-exp(-^M) 



\si-Sj\> eT 1 /' 3 . 



(4.34) 

Therefore, for any given k' > 0, by taking e > small enough but fixed 
and then taking To large enough, there is a constant C > such that for all 
T>T , 

'T\ 1/3 Z(j). 



xj Z(i) 



—Vijiu^Uj) 



< C • exp(— n'\si — s 



3\)i 



(4.35) 



and hence 

'T\ 1/3 A(j) 
x) Aii) 



Vij(ui,Uj) 



< Cexp(— K ! \si — Sj\ + TjSj — TiSi). (4.36) 



Finally, given any k > 0, by take k! > n + max{|ri|, • • • , \r m \}. Then 



when Si — Sj > 0, then 



TjSj - TiSi + (K - K)\Si - Sj\ = -(n - Tj)Si - Tj(Si - Sj) + (k - K?)\Si - Sj 



< -{n - Tjfa < -~{n - Tj){ Si + Sj) 



(4.37) 
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since T; — Tj > 0. Similarly, when Sj — Sj < 0, then 

T j S j T i S i "I" ( K K )\ S i S j\ = T j) S j s j) + ( K ^ )\ S i S j I 

(4.38) 

This implies (Q4"]l . □ 
We now prove an asymptotic result for X. 

Lemma 4.4. Consider the scaling fl^.ffp . T/ien, 
'T\ 1/3 A(j) 



-K itj (u h Uj) 



(4.39) 



1/3 N 



/■oo 

= / c?A Ai(A + Si + if) Ai(A + Sj - + r?)e- A ( T '-^ + 0(T" 

uniformly for Si, Sj in a bounded set. 

Proof of Lemma 14.41 By the definition ( 12. 14ft . 

T\ 1/3 ~ . . /T\ 1/3 -1 / , / , e-^^(w) 1 

(4.40) 

The steepest-descent analysis of integrals very similar to this one with the 
same scaling (14.21) has been performed repeatedly in various places (see for 
example, [21IH1[IIJ[]]3[IH])- The only difference here is that we have a double 
integral and we have to make sure that the two path do not touch. However 
this can be easily handled by locally modifying the steep(est)-decent contours 
near the critical point. Except for this modification, the analysis of our case 
is similar to those in the literature. Nevertheless we provide the proof for 
the completeness of the paper, and also since the analysis of this Lemma 
and Lemma 14.51 is a prototype of all the other Lemmas in this section below 
(except for Lemma I4.9p . 

The first step is to determine a steep descent path for <pi(w)e~ WUi and 
(f)j(w)~ 1 e ZUj . Let us write 

(j>i{w)e- WUi = exp {Th (w) + T 2 ' z h hl {w) + T^h^w)) , (4.41) 
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where 

h (w) = -w + p 2 ln(l/2 + w)-{l- p) 2 ln(l/2 - w), 
h lti (w) = r lT ±— [(1 - 2p)w - x ln(l/4 - w 2 )] , (4.42) 

Note that ho(w) is independent of i. For the steep descent analysis we need to 
determine a steep descent path for ho. The steep descent path will always be 
taken symmetric with respect to complex conjugation. Thus we can restrict 
the discussion below to the part lying in the upper half plane. We have 

h n (w ) — — 1 + — 7- h — , hJw)- 



0V ; 1/2 + w 1/2-tu' 0V ; (1/2 + w) 2 (1/2 -w) 2 ' 

(4.43) 

which both vanishes at the critical point u> c = a = p— 1/2, and h' "(w c ) = 2/x- 
Let 7^ = {w = i - ae ; V e [vr/2, 3vr/2)} with a > 0. Then, for 
6 1 G [7r/2,7r), Re(/io(w)) is strictly decreasing in 6 1 . Indeed, 



A. 
d§ 



Re(h (l - ae w )) = -a sm(9) (l - r^j^ ) < 0, (4-44) 



since the last parenthesis is strictly positive: for 6 G [n/2,7c), Re(w) > 1/2 
and \w + 1/2| > 1, while p < 1. In the neighborhood of the critical point we 
consider a second path, = {w = p — 1/2 + e~ 17r / 3 t, t G [0, 2(1 — p)]}. Along 
that path, we have 

l R e (ftoH ) = - g^rt + P^+f) , (4 . 45) 
at 2 1 1/2 — tw| 2 |l/2 + 

The second degree term, 2p(l — p) + (1 — 2p)t + 1 2 , is strictly positive for all 
p G (0, 1) and t G [0, 2(1 — p)\. Therefore also 71 is a steep descent path, and 
close to the critical point will be steepest descent. 

Now we can define the steep descent path used in the analysis. Let 
7 = 71 U 7^/3( 1 _ p ) U 71. In a similar way, one obtains a steep descent path for 
— ho(z), namely, T is the path obtained by rotation around the origin of the 
steep descent path 7 but with 1 — p instead of p. Finally, as we shall discuss 
below, any local modification of the contours in a region of order T -1 / 3 around 
the critical point is allowed, see Figure |4TT| so to have \z — w I > eT" 1 / 3 for 
a fixed e > 0. 
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Figure 4.1: The modified contours 7 and T are both anticlockwise oriented. 
The modification of T close to the critical point a is of order T~ 1//3 . 



The path 7 and V are steep descent path for ho(w) and —h (z). Therefore, 
for any given small 5 > 0, the contribution to the double integral ( I4.40P 
coming from 7 xT\{|;z—a| < 5, \w— a\ < 8} is only of order Z(i)/Z(j)0(e~ flT ) 
for some /i = fi(5) > (with /x ~ 5 3 for small 5), which is smaller than 
Z{i)/Z(])0{T- 1 ^). Note that Z(i)/Z(j) is the value of the integrand at the 
critical point. 

Next we analyze the contribution coming from a 5-neighborhood of the 
critical point w c = a. There we can use Taylor series expansion of ho, hi ti , h 2 ^, 
which are given by 



ho(w) = h (w c ) + -x \w - w c ) 3 + 0((w - w c ) 4 ), 
hi,i(w) = h lti (w c ) + n X ' 2/3 ( w ~ Wc? + 0(n(w - w c ) 3 ), 



(4.46) 



Therefore, the main contribution to (14.401) is given by 



Z(i) (T 



1/3 



-1 




e T(w-W c ) 3 /3 X +T 2 rA Ti{w _ Wc )2/^/3_ T l/-A Si{w _ Wc)/x l/3 



z(j) \x 



(27Ti) 2 



dw 



e T(z-w c )y3x+T 2 / 3 T 3 (z- Wc ) 2 /x 2 ^ 3 -T 1 /3 Sj (z-w c )/x 1 / 3 
^(T^-^V^ 2 / 3 ^-^) 3 ) 1 




33 



where 75 and T$ are the pieces of 7 and T which lies in a (^-neighborhood of the 
critical point w c = a. More precisely, setting Hi := h + T~ 1 / 3 hi i + T~ 2 / 3 h2,i 
one has that 

THi(w c +t)- TH t (w c ) - Thf{w c )^ - - T 1/3 -^t 

<T sup |4 4) HI|t 4 |/4! + T 2 / 3 sup I hf\ (w) ||t 3 1/3!. 



B(w c ,\t\) B(w c ,\t\) 

Assume that < 77 < | min{p, 1 — p}, then 



(4.48) 



and 



sup |/^>)l < 6.2 4 (p"^ + (1 - p)^) (4.49) 

B(w c ,ri) 

sup l/i^HI < N^4V 3 + (1 - p)- 3 ). (4.50) 

Thus, it is an easy computation to check that one can find r\ > small 
enough so that 

r? sup \ht\w)\<hf{w c )%- and sup |^Wh 3 < |^> c )|^. 

(4.51) 

Assume that < 5 < rj. Then it is easy to show as in (14.22ft . Lemma 14.21 
that 

z(i) fr\ 1/3 -1 f dz r dweTHiiw) _ THj{z) _ 1 



Z(j) \xj (2vri) 2 Jr 6 J ys ' w-z 



Z(j) Vx7 (2vri) 2 ~~ ^ '*' T^^+TVs Z^f^ -TVafi^ w-z 
< C(T~ 1/3 ). 

(4.52) 

It thus remains us only to determine the asymptotic of (14.521) . We make 
the change of variables W = (w — w c )(T/xY^ 3 and Z = (z — w c )(T /x) 1 ^ 
and obtain 

7(i\ -1 /•e 2 ' ri /3 < 5(T/ x )V3 r e-^/38(T/ x y/3 p W 3 /3+TiW 2 ~ Si W 1 

' dZ dW * 73/ ^ 72 _ 7 J - (4.53) 



(27ri) 2 J e -w s{T/x) i/3 7 e .i/3 5(T/x)1 /3 e^ 3 /3+^^- Sj z ^ _ z 
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where the integration paths do not touch. We can now extend the integration 
path to infinity, since 

7(i\ _1 /•e 2 ' ri / 3 oo r e-^/ 3 oo w 3 /3+n W 2 - Si W i 

1 dZ dW 



Z{j) (2vri) 2 y e - 2 . i/3oo Jewaoc eZys+TjZi-s.z w _ z 

Z(i\ _1 re 2 ^/ 3 8(T/ x ) 1/s re-^ 3 8(T/ x ) 1/3 e W 3 /3+ n W 2 -s,W ]_ 
W)W? Je-^S { T,J 3 Je^ 3 8(T/ x) ^ ^ W - Z 

< Ce- 53T/ ^. 

(4.54) 

The conjugation by A(J)/A(i) and the identity (|A.2[) end the proof. 



□ 

Now we give a bound which holds uniformly for Sj, s 3 - bounded from below. 
Let So G R be given. 

Lemma 4.5. Consider the scaling RJfy . Then, for any given k > 0, there 
exists a T > /arge enough such that for all T >T and Sj, s 3 - > s 



^ Mi) „ <Ce-<"+">. (4.55) 

The constant C depends on k, Tj, t j only. 

Proof of Lemma 14.51 The upper bound of the integrals similar to (I4.40p 
has also been obtained in various places (see for example, [2|l6tlTT |[T7l[T9] ). 
The analysis in our case is similar to those in the literature. However, again 
we provide the proof for the completeness of the paper; also the analysis in 
this proof is going to be used and adapted in all the Lemmas below in this 
section (except for Lemma [4. 9p . 
First of all, we can rewrite 



Ki,j(ui, u 



•3) 



— , i d\[ — (b dw6i(w)e 1 ?^ (h dz 

(4.56) 

Set uq = T — Ti 2( ' 1 ^^x — (which corresponds to Sj = 0). Then 

1/3 



-) ^Kijfauj) = [ dXEtisi + X)E 2 { Sj + A), (4.57) 
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where 





r 








V2vri 




r 








V27ri 



' Fl/2 1/3 ( 4 - 58 ) 



We now show the exponential decay of E\(si) for large positive Sj. 
Lemma [4.41 indeed implies the result when sq < Sj < 0. It is thus enough to 
consider the case where Sj > 0. The analysis of E 2 (sj) is made in exactly the 
same way, up to a rotation around the origin and exchange of p with 1 — p. 

To this aim we modify the contour given on Figure |4~T1 as follows, defining 
a new contour 7'. Call 7'+ the part of the contour lying in the upper half- 
plane {Imz > 0}. Let k > be a fixed constant and 

e(i) := 2\n\ + 2k. (4.59) 

Then we set 

7 '+ = { a + feW*, t > £ (i)(T/ X )- 1/3 } U {a + £(i)(T/ X )- 1/3 e ie , < 9 < vr/3}. 

(4.60) 

The contour 7' is then completed by adjoining the conjugate of 7 /+ . This 
modification of the contour 7 has no impact on the saddle point analysis 
made in the proof of Lemma 14.41 Indeed, for T large enough 

\Th (a) + T^h^a) - Th (a + e(t)(T/ X )- 1/3 e w ) 

- T 2 / 3 h hi (a + £(i)(T/ X )- 1/3 e w )| < C(i), (4.61) 

where C(i) = 2e(i) 2 (\Ti\ + e(i)) is a uniformly bounded constant as 7$ is 
chosen in a compact interval of R. 
Along 7' it holds, 

Re(w -a)> ^(T/x)- 1 / 3 > (« + MXr/*)" 1 /*. (4.62) 



Thus, along 7' 

' < e- KSl . (4.63) 



e -(v>-a)<>i{T/x) 1/a +Ti8i 
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(a) (b) (c) 

Figure 4.2: Contours used in the following Lemmas. The paths have local 
modifications close to the critical point a that are only of order T -1 / 3 . 



Using the saddle point argument used in the proof of Lemma [4.41 and (|4.61|) . 
we easily deduce that there exists C > independent of such that for T 
large enough 

< Ce~ KSi . (4.64) 



Thus, 



POO 

|(!457D| < / d\\Ei( Si + X)\\E 2 ( Sj + A) | < Ce-«(' i+ 'A (4.65) 
Jo 



for some other constant C . 



□ 



The saddle point analysis made for K in the preceding Lemmas will now 
be used, up to minor modifications, to consider the asymptotics of the re- 
maining terms. 

Lemma 4.6. Uniformly for s\ in a bounded interval one has that 



lim 

T-s-oo 



1/3 



Ra, 



_ a = Sl + e-H— 



1 dx dy A\{rl + Si + x + y)e' 
Jo Jo 



n (x+y) 



(4.66) 



Moreover, for any given k > 0, there exists a T > large enough, such that 
for all T > To it holds 



T 



\ 1/3 



XJ 



j Ra-a — S\ 
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< Ce 



-KSl 



(4.67) 



for any s > Sq. The constant C does not depend on T or s\. 
Proof of Lemma 14.61 We have 



TV 1/3 fT\~ 1/z 1 / , (hAaje 

and also 

,,1 M = exp [-(TV*) + T^i,!^) + T^h^z))} (4.69) 

with hkS defined in (14. 46 p . As steep descent path we use almost the same 
contour T as in Figure I4TT1 up to T -1 / 3 -local deformation so that it passes to 
the right of a, see Figure fl~2l (a). Then, using the same asymptotic argument 
as in the proof of Lemma 14.41 we deduce that 

&U) R ~ = ™J I , dZ — * — • (4 ' 70) 

where with r ) we mean the path from e~ 2n ^ 3 oo to e 27T ^ 3 oo and passing on 
the right of 0. We can take it to its left up to adding a residue term, which 
is just S\. Then the identity ( ]A.4j) gives the first formula of the Lemma. 
For the exponential decay, one first computes the residue at z = a, yielding 
that 

R a,~a = Si + - — f dz * 1 ) ~ rj. 4.71 

The contour now lies to the left of a, which ensures using the same argument 
as in Lemma 14 .5} the exponential decay for large positive s±. 

□ 

Now we we turn to the asymptotics of the function g defined in 02. 1 2j) 
and the different terms of the function F defined in (13. 7J) . 

Lemma 4.7. Consider the scaling Q3ty . Then 

lim }fl 2 \ Ui ) = - dy Mtf + 8i + y)e™ (4.72) 

T^oo A(l) J 
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uniformly for s, in a bounded set. Moreover, for any given k > 0, there exists 
a T > large enough, such that for all T >T it holds 



1 



A(i] 



(1/2), 



U: 



< Ce' 



(4.73) 



for some constant C independent ofT and uniform for Sj bounded from below. 
Proof of Lemma 14.71 We have 

e Th (w)+T 2 / 3 h 1>1 ,(w)+T 1 / 3 h24w) 



(1/2), 



'It; 



1 

27Ti 



r i/2 



w — a 



(4.74) 



with the hkS defined in ( I4.46p . As a steep descent path we use 7 in Figure H~T1 
with a local modification on the T -1 / 3 scale so that it passes on the right of 
a, see Figure I4T21 (b). Then, the same asymptotic argument of the proofs of 
Lemma 14.41 and Lemma 14.51 still holds and gives us the claimed result. We 
also use the identity (IA.5j) to express the result in terms of Airy functions. □ 

Lemma 4.8. Consider the scaling \4-%fy . Then 



lim 



T->oo A(i) J U1 



dyK hl (u i ,y)fi a) (y) 



e-a^f-n-i / d\ dye 



-A(n-Tj) -nv 



e~ TlV Ai(r/ + Si + A) Ai(r 1 2 + Sl + A + y) 

(4.75) 

uniformly for Sx,Sj in a bounded set. Moreover, for any given k > 0, there 
exists a Tq > large enough, such that for all T >Tq it holds 



1 



J U1 



dyK i)1 (u i ,y)fi a) (y) 



< Q g-«(»i+«i) 



(4.76) 



for some constant C independent of T and uniform for Si,s\ bounded from 
below. 

Proof of Lemma 14.81 Notice that f[ a \y) = 4>i{a)e~ ay . Then, choosing the 
integration path for z in ( I2.14p so that Ke(z) < a, we can explicitly integrate 
the variable y and get 



dyK iil {u h y)f[ a) {y) 



(27ri) 



dz 



r_i 



/2 



dw 



r i/2 



e zux ~ wu ^i[w) 4> 1 (a)e- aui 
(w — z)(f>i(z) a — z 

(4.77) 
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where the integration path for z does not include a and does not cross the in- 
tegration path for w. The rest follows along the same lines of Lemma |4~41 and 
Lemma T4.51 Finally one uses (1A. 6|) to rewrite the double integral expression 



in terms of Airy functions. □ 
Lemma 4.9. Consider the scaling fiJJty and let % > 2. Then 



J im tt\ / d y v iA u i>y)fl (y) = —r== , , / d v e Mn ~ Tl) - 

(4.78) 

uniformly for Sj — s\ in a bounded set. Moreover, for any given k > 0, there 
exists a T > large enough, such that for all T >T it holds 



(4.79) 

for some constant C independent ofT and uniform in s$, Si. 

Proof of Lemma 14.91 Using /] (y) = 4>i(a)e~ ay and modifying the inte- 
gration path in ( 12.141) to iK. + a for any a e (a, 1/2), we can integrate out 
the y variable with the result 

dyV iyX {ui,y)f[ a \y) 



Zi) 



1 1 I dze ^tMM«)e-^ 



z (i) 2vri J iR+a ' 0i (z) z-a 



2vri J iR+a " 0i(*)z-a 

For the saddle point analysis we then deform back our integration path to 
a + iM with only a local deviation order T~ 1//3 around z = a to make it passing 
on the right of a. Specifically, for a given small e > 0, we use the integration 
path 

7 = {d + ix,iGt\ (sT- 1/3 , eT- 1/3 )} U {a + eT" 1/3 e ie , 9 e [-tt/2, tt/2]}, 

(4.81) 

see Figure H~2l (c). The asymptotic analysis for - Si in a bounded set is 
the same as in Lemma [4.21 except that the integral has an extra 1/Y in the 
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denominator. Explicitly, we get 



1 /•«! 1 /" p-^ 2 (r i -ri)-iy(s i -si) 

lim _ / dyV t , 1 (u l ,y)ft\y) = -- dY 



for any e > 0. 
The function 



iy 

(4.82) 



1 f -Y'in-r^+iYs 

v(s) :=— dY (4.83) 

2n ./to i~ ir 



rf w(s) = 1 exp [—77-^- r ) (4.84) 



satisfies 

d~s V ^ s)= ^(n-n)^ V 4(ri-n). 

and the boundary condition 

lim = 0. (4.85) 

s— >— oo 

Therefore, t> (s) is a Gaussian distribution function, so that our limiting func- 
tion is 

1 f Ul ( a ) 1 rsiSi 

lim —— \ dy Vn(ui, y) U a) (y) = — = / dxe 4 K-n). 

(4.86) 

In order to prove (14.791) . we adapt the proof of Lemma I4T31 First, suppose 
that Si — Si > 1. For this case, we take the contour 7 = {z c + it\t G R} as 
in the proof of Lemma H~3"l Note that dist(a, 7) = dist(a, z c ) > —, — x ^ 1/3 . 
Hence and the term in (I4.80P is bounded by 0(T 1 ^) and an analysis 
as in the proof of Lemma 14.31 as in the proof (noting the presence of T 1 / 3 
in (|%2U ) implies (H7T9]) . 

Secondly, suppose that Sj — s\ < —1. We take the same contour 
7 = {z c + it\t G K}. In this case, the residue at the simple pole a should 
be taken into account. Except for this term, the analysis on the contour 7 is 
the same and we obtain (14. 79 p . 

Finally, when — 1 < Sj — s\ < 1, we take the contour 7 = 71 U 72 where 
71 is the part of the contour {z c + it\t G M} which lies outside the circle of 
radius r := ^ _ r ^ Tl/3 centered at a, and 72 is the part of the circle of radius 
r centered at a whose real part is at least z c . This contour is same as the 
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straight line {z c + it\t G R}, except that it goes around a on the right in 
a C(T _1 / 3 )-neighborgood of a. On 7 1; as dist(a, 71) > r = C(T~ 1//3 ), and 
hence we obtain, as in the proof of Lemma I4.3[ 

Z(l) 1 



Z i) 2ni 



dze 



zim-m) <t>i{ z ) 1 
0i(z) 2 — a 



< e - K l s ;- s il 



(4.87) 



On the other hand, on 72, observe that dist(a,72) = 0(T~ 1 ^) and the arc 
length of 72 is also 0{T~ 1 ^). By using the estimate (I4.18P and (I4.19p . which 
holds on 72 if T is large enough, we obtain that T 2 ^ 3 g(z) is bounded for 
z G 72 (where g(z) is in (14.7[) ). Hence (recall the first equation of (I4.26P ) 



Z(l) 1 



Z(i) 2ni 



dze z( - Ul ~ Ui) - 



72 



\\z) z — a 



Tyz\g{z)-g{a)\ 



72 



2 — a 



This, together with (Q7|l implies (STSJ. 

The last term to be computed is the asymptotic of g. 

Lemma 4.10. Consider the scaling (fjJ?D. Then 



0(1). 
(4.8<' 



□ 



lim A(j)gj(uj) 

1 —too 



POO 

- dx Ai(r ? 2 + Sj + x)e~ T ^ 
Jo 



(4.89) 



uniformly for s j in a bounded set. Moreover, for any given k > 0, there exists 
a Tq > large enough, such that for all T > Tq it holds 



\A(j) gj ( Uj )\ < e a^+^ + Ce-^ 



(4.90) 



for some constant C independent of T and uniform for Sj bounded from 
below. 



Proof of Lemma 14.101 We have 



9jW> 



1 

27Ti 



dz- 



(Th Q {z)+T 2 ^h^ ( Z )+T 1 ' 3 h 2 tj (z)) 



r_i/2 



z — a 



(4.91) 



with the /ifc's defined in (I4.46p . As steep descent path we use T in Figure H~T1 
with a local modification on the T _1//3 scale so that it passes on the right 
of a, see Figure 14.21 (a). Then, the asymptotic argument of the proofs of 
Lemma [4.41 and Lemma [4.51 still holds and gives us the claimed result. The 
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only small difference is that for the bound the integration path for large s/s 
has to be chosen to pass on the left of a, which is fine up to adding the 
residue at a. The factor e 2T ^ 3+TjSj comes from the multiplication by A(j). 
To rewrite the result in terms of Airy functions we used flA.3j) . □ 



5 Proof of Main Theorems 

In this section we prove Theorem 11.21 Recall (14. ip 

f]{G(x(r k ),y(r k )) <£(r k ,s k )}\ =J2^~ k ' det ( 1 ~ P ^ P t 

k=l / k=l 

(5.1) 

where 

Qt ■= ( |J {Ra,-a - ((1 - PJCPt)- 1 P l F,P t g)). (5.2) 

and a = p — 1/2. We take the limit T — > oo with the scaling (II. 5p . 

Note that the scaling of (I1.5P is related to the scaling (14 .2p by the identities 
\x{ji) + yiji) — 2t\ < 2 and \xiji) — y(rj) — 2t^\ < 2. The difference (at most) 
±2 does not contribute to the asymptotics and only makes the notations 
complicated. For this reason, we will ignore this difference ±2 in the following 
presentation. 

Proposition 5.1. Fix m G N ; and real numbers T\ < t 2 < . . . < T m , 

s±, . . . , s m G M. Let x(r k ), y{r k ) and £(r k , s k ) be defined as in (II. 51) . Then it 
holds 

lim det (l - P e K c ° ni P e ) = det (t - P s K Ai P s ) 

T^oo \ /L 2 ({l,...,m}xK) V /i2({l,...,m}xM) 

(5.3) 

where P t (k,x) = l[ x> e( Tk ,s k )] and P s (k,x) = t[ x>Sk }- 

Proof of Proposition 15.11 Proposition 14.11 implies that for Sj, Sj in a 

bounded interval / 



lim ( — ) ^^K^(ui,Uj) = [K Ai \ i;j (s i ,s j ), (5.4) 



"T\ 1/3 A(j) . 

and that the convergence is uniform on /. Now Lemmas 14.31 and 14.51 
then imply the convergence in trace class norm of the rescaled kernels 
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(T/ X ) 1/3 ^K ltJ (u u u,) and {T /xf^^V^u,). This completes the 
proof. □ 

Proposition 5.2. Fix m G N real numbers T\ < t 2 < . . . < r m , 

Si, . . . , s m G R. Lei x(rfc), w(Tfc) and £(7%, s^) 6e defined as in (jl.5p . T/ien 

lim ft = g m {r,s) (5.5) 

T— 5>oo 

where g m (r,s) is defined in (11.1 ip . 

Proof of Proposition 15.21 Consider the term ((1 — P t K Pe)^ 1 PeF, P e g) . 
We first scale and also conjugate by some operators. Define the function 

:=T — Ti 1-2 ^ T / + a; (5.6) 

Set /5 = 1 + max \t{\, and define the multiplication operator 

i=l,...,m 

W(z,a:) :=^(i)e"^ x . (5.7) 

Set 

Ll(x,x') := (£) 7 W^x^jiu^^ix'^W-^x'), (5.8) 

and 

*f (x) := W(i, x)g i {u i {x)), $J(x) := W 1 ^, x)F i {u i {x)), (5.9) 

where one recalls the definitions of Fi in (13. 7p and g>j in ( I2.12p . Here the 
superscript T indicates the dependence on T. 

By Lemmas 14.71 14.81 14.91 and 14.101 for any given k > 0, there exists a 
T > large enough, such that for all T > T , tyf and $f are £ 2 ([sj, oo)) for 
any given Sj. Then for T large enough, 

/ \ 1/3 

1 1 J ((1 - PeKP^Pt^Ptg) = ((1 - P a L T P.)~ 1 W T ,-P.* T )- (5-10) 
Let Q be the multiplication operator Q(i,x) = e~^ x and set 

p=(l-P s Q^AiQ~ 1 i , »)~ 1 , $ = Q*, $ = Q _1 $ (5.11) 
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where the functions are defined in (11.81) . Observe that the functions 
e L 2 (( Si , oo)). Since 

((1 - l' s f\\]l\ ) I'Ak P s ^> = (pP s $, P s $>, (5.12) 



it is enough to prove \im.T^oo{p T P s & T , P S ^ T ) = (pP s &,Ps^) where 
P 1 

Clearly, 



•T-'I-P.LTPS- 1 



(p'P.Q^P.V^-tfP&P.V) 



< ||p T -p||||P s $ T ||||P^ T | 



+ ||p|| ||P s (^ t - *)|| ||P S $ T || + ||p|| ||P.*|| ||P s ($ t - $)||. (5.13) 

First consider \\p T — p\\. As the operator-norm is bounded by the Hilbert- 
Schmidt norm II ■ II 2, 



\P S L T P S - P s QK M Q~ l Psf < \\P S L T P S - P s QK M Q- l P s \ 

m POO POO 

Y, dx dx ' \ L U X ' x ') - e-^ x - x,) [K Ai ] ltJ (x, x')\ 

J Si J S-i 



By taking k = + 1 in Lemmas 14.51 and I4.3[ and using the dominated 
convergence theorem and Proposition I4.1[ the limit as T — > 00 of f !5.14j) 
becomes 0. By Lemma IB. H the operator 1 — P s Ka_\P s is invertible and it 
follows that \\p T — p|| — >■ as T —> 00. 

Now consider the term P S ^ T . From (I4.90p of Lemma 14.101 we have 
|^J(^)| < Cie~^~^ x + C 2 e~ l - l3+ti)x for x > Sj for a fixed Sj. The conjuga- 
tion e~ l3x is introduced to make this function decay exponentially as x — > 00. 
Hence ||P S , I /T || is bounded uniformly for large enough T. Also noting that 
the right-hand-side of (I4.89P is precisely \F,(s,,), it follows from dominated 
convergence theorem that ||P S (\1/ T — *ff)\\ — > as T — > 00. 

It follows from Lemmas S21 SSI and S3] that \$j(x)\ < Ce~^ x for 
x > Sj for a fixed Sj. Hence ||P S $ T || is uniformly bounded for large enough 
T if we set for example k = (3 + 1 > 0. Also it follows that ||P j ($ T — $) || -> 
as T 00. 

Hence we have shown that (15.1 Oft converges to (15.121) . The remaining 
term in Qi is (x/P) 1//3 P Q .-a- But Lemma [4.61 shows that this converges to 1Z 
(after the changes of variables x 1— > x — s\). This completes the proof. □ 
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The above two Propositions proves that Qt ■ det (l - P/iT° nj P £ ) con- 
verges to g m (r, s) det (^1 — PsKaiPs^ as T — > oo. It remains to show 
that the derivatives with respect to s^'s also converge. Let us first con- 
sider (x/T) 1 / 3 ((l - PeKP^P^Peg) of Q t . As in floTTUj) . this equals 
((l-P a L T P a )- 1 P a $ T , P S ^ T ). Note that the dependence on Sk is only through 
the projection operator P s . Hence by simple translations, 

((1 - P s L T P s )- l P s ^ T , P^ T ) = ((1 - L T )~ 1 <fr T , * T ) (5.15) 

where L T ((i, x), (j, y)) := L T ((z, x+sj, (J, y+Sj)), $ T ((z, z)) = $ T ((i, z+s f )) 
and ^ T ((i, x)) = ^ T ((i, x + Si)), and all the (real) inner product, the operator 
and the functions are all defined on L ({1, . . . , m} x R+). Then 

_ L T ) _1 * T , * T ) = ((1 - L T )-\^--L T )(l - L T )- 1 * r J * T > 

+ ((1 - L T )- 1 (A$T ) + ((1 _ jTyltfT JLyTy 

OS k OSk 

Hence we need to control the asymptotics of the derivatives of the kernels. 
But since all the asymptotic bounds for large x in the Lemmas in SectionSJare 
all exponential, while the derivatives of the kernels yields only polynomial 
terms, we can check that the dominated convergence theorem still applies 
and obtain the convergence of (15.161) to the corresponding derivatives of the 
limit. We omit the detail. The other terms are also similar. This completes 
the proof of Theorem 11.21 

To prove Theorems II .541X771 we use the following slow-decorrelation result, 
which is an extension of Proposition 8 of 



Lemma 5.3. Let A = (c\T, C2T) for some C\,C2 > 0. Let then 
B = A + r((l — p) 2 , p 2 ) with r ~ T v with < v < 1. Then, for any 
13 e (i//3, 1/3), it holds 

lim P (\G{B) - G(A) -r\< T?) = 1. (5.17) 

Proof of Lemma 15.31 First of all, it is easy to check by Proposition 2.2 
of [T7] that the difference between the model (II -ip and the last passage per- 
colation model corresponding exactly to the stationary TASEP becomes ir- 
relevant in the order as T — > 00. Denote Pta the measure for the last passage 
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percolation corresponding to the stationary TASEP. Due to the stationarity 
in space and time of TASEP, we observe that G(B) — G(A) and G(B — A) 
have the same distribution. Therefore, 

F TA (\G(B) - G(A) -r\< = ¥ TA (\G(B - A)—r\< T?) 

= P TA (|G(r(l-p) 2 ,rp 2 )-r| <T^). 

(5-18) 

But since the distribution of (G(r(l — p) 2 ,rp 2 ) — r)/(r/xY^ 3 converges to 
F pi], and r^/T 13 = 0{T P ^ P ) as T -»■ oo, we find that f l5T8|) 
converges to 1 as T — > oo. This completes the proof. □ 

Proof of Theorem II. 5L The projection of (x(r,9),y(T,8)) on the line 
{x + y — (1 — 2x)T} along the critical direction ((1 — p) 2 ,p 2 ) is (x(r),y(r)) 
with x, y defined in (11. 5p . We have 

x(t) = x(r, 6) - r(l - p) 2 , y(r) = y(r, 6) - rp 2 (5.19) 

with r = r(6) = 9T U . In particular £(r, 0, s) = £(t, 8, s) — r(8). Then, 



P (f]{G(x(T k ,9 k ),y(T k ,9 k )) <l(r k ,9 k ,s k )} 



\k=l 



(5.20) 



= P f]{G(x(T k ,0),y(T k ,0)) <£(r k ,0,s k ) + E k } 

\k=l 

with 

E k := G(x(r k , 0),y(r k , 0)) + r(9 k ) - G(x(r k , 9 k ), y(r k , 8 k )). (5.21) 
By Lemma [5.31 we have S k = o(T 1 / 3 ) so that 

/ m 



lim(E20])= lim P( C]{G{x{r k ,0),y{T k ,0)) < £{r k ,0,s k )} , (5.22) 

\k=l 



see the proof of Theorem 1 of [15] for detailed steps. This is by Theorem 11.21 
the desired result. □ 

Proof of Theorem II. 6L The goal is the express our setting into the one of 
Theorem 11.51 Denote p(r) := q(r) + n(r), which is given by 

p(r) = (1 - P ) 2 T + r2(l - p) X 1/3 T 2 / 3 - (1 - p)sT^/ x ^ (5.23) 
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and recall 

n(r) = p 2 T - t2 PX 1/3 T 2 / 3 . (5.24) 

(Since n(r) and q(r) are integers, the above formulas are are exact up to the 
error of size 2 at most: since this difference does not affect the asymptotics 
but only complicates the formulas, we drop this difference in the following 
presentation.) By (11.161) we have 

P (f]{K n(Tk) (T) > g(r fe )}j = P (f]{G(p(T k ),n(r k )) < T}j . (5.25) 

Denote by (X(t),Y(t)) the projection of (p(r),n(r)) on the line 
{x + y = (1 — 2x)T} along the critical direction ((1 — p) 2 , p 2 ). We get 

X(t) =p(r) -r(l -p) 2 , Y(t) = n(r) - rp 2 (5.26) 

where 

2(1 - 2p) y 1 / 3 T 2 / 3 1-p fml . 1/3 r _ _, 

r = r r,s = r ; « s- — (T/x) 7 ■ 5.27 

1 - 2x 1-2% 

Let us further denote 

T ' = T ~ ^n rm - (5 - 28) 

Then, 

T = 1(t s ,s) +r(r,s) (5.29) 
where £(r, s) is defined in (11. 5ft . Then, replacing ( 15.271) into ( 15.261) we get 

X(r) = (1 - pfT + r 2x4/3r2/3 - s -^It 1 / 3 = a;(r s ), 

X " 2 ^ 1 ~ 2X (5 30) 

2y 4/3 T 2/3 2/3 ^' OU ^ 

F t =pT- t— — + s - —T /A = y{r , 

1 - 2x 1 - 2x 

where x(r),y(r) are defined in (jl.5p . 

With these notations we can rewrite ( I5.25P as follows, 



(JEIBJ = P f]{G(x(r s k ) + r k (l - p) 2 , y{r s k ) + r fc p 2 ) < £(r fc s , s k ) + r k } 



\k=l 



(5.31) 
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where := r(rfc,Sfc) ~ 0(T 2 ^ 3 ). Then, by Theorem 11.51 the result follows. 

□ 

Proof of Theorem 11.71 The proof of Theorem 11.71 is very similar to the 
one of Theorem 11.61 The scaling (jl.22p corresponds, in terms of x, y, to the 
scaling 



x = (1 - p) 2 T + 2r(l - p) X 1/3 T 2/3 - sx^T 1 / 3 , 
y = p 2 T- 2rp X 1/3 T 2/3 - s X 2/3 T 1 ^. 



(5.32) 



The projections (X(t),Y(t)) of (x,y) on the line {x + y = (1 — 2%)T} 
along the critical direction ((1 — p) 2 ,p 2 ) are given by X(t) = x — r(l — p) 2 , 
^( r ) = V ~ r P 2 - Explicitly we find 

, * 2(1 - 2p) X 1/3 T 2 / 3 2 X 2/3 T 1 /3 

r(r ' s) = r T^c S T^2F' 

X(,)^ 2 T + ,g^ + g ( 1 - 2 ^ 2/3TV3 , (5.33) 
K ' P 1-2* 1 - 2 X 

2 X 4 / 3 T 2 / 3 (1 - 2p) X 2 / 3 rV 3 

it = pi — r s . 

1 ' 1-2* 1-2* 

r'^M^l (5,4) 



Setting 



we get 



T = £(r s ,s)+r(r,s), X(t)=x(t s ), Y{r) = y{r s ). (5.35) 



Then, setting r k := r(r fc , s fc ) we get (I5.3ip and then by Theorem [L5] the result 
follows. □ 



A Some Airy function identities 

In the asymptotic analysis we get two basic integral expressions which can 
be rewritten in terms of Airy functions and exponential, namely 



-1 
2vri 
1 
2vri 



/ dWe WV3 +b Wi-cW = M{b 2 + c)e |6 3 +te 

' x (A 1) 

dZe -Z*/ 3 - b Z* + c Z = Ai(6 2 + c)e - - 



'oo 

; 2^/ 3oo 
-27ri/3 c 



49 



Starting from these two formulas we state some identities. 
Lemma A.l. 

„ e 27ri/3 00 n e - ni / 3 oo W 3 /3+b 1 W 2 -c 1 W i 

( A ) 7^2 dZ dW 



(27Ti) 2 J e -^ /3oo 7 eW3oo e Zy3 + b 2 Z*-c 2 Z W _ Z 

e |fcf+6ici /-oo 



ff>3+b 2 c 2 



e 3 



/•oo 

/ dXe-^ b2 - bl ^ ki(b\ + Cl + A)Ai(6 2 , + c 2 + A), (A.2) 



i r e 2ni/3 oo , 
1 / dZe _ z s /3 _ bz * +cZ ^ 

poo 

/ dxAi(i 2 + c + i)e" k , (A.3) 
Jo 



(B) ,_. , 

Z7T1 ./ e -2vri/3 00) ^ o/ Q 



e -§b 3 -(. 



-, r e 2ni / 3 oo -, 
/7T1 J e _2^i/ 3o05 ^ /o Zj 



nOO POO 



poo poo 

dx dy Ai(b 2 + c + x + y)e- b{x+y) , (A.4) 
Jo Jo 



i re "'i-oo i 

(D) — / dWe z*/3 +b w 2 -cwJ_ 

/•oo 

,c / diAi(l) 2 + c + i)e k , (A.5) 



' e 7Tl /' A oo, right of 

2 

= es 



(E) — — / / 



/»oo /»oo 

/ dX / ^e- A(62 - 6l) Ai(6 2 + Cl + A)Ai^ + c 2 + A)e- 62X , 
Jo ./o 

(A.6) 



(F)For6 2 <6i, / d\ e^ 2 ^ ki(b\ + c x + A) M{b\ + c 2 + A) 



!> I ^ + 1(4 - bl) + b 2 c 2 - b lCl ) . (A.7) 



v/4tt(6i - 6 2 ) V 4 ( fe i - fo 2) 3 
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Proof of Lemma IA.ll (A) Since we can choose the paths for W and Z 
such that Re(VT — Z) > 0, using 



1 



W-Z 



d\e~^ w ^ (Ai 



o 



and (TO) we get flA~2l) . 

(B) Let /(c) := l.h.s.( 1A.3"l) . Then differentiating and using 11A.1I) we 
get f'(c) = — Ai(6 2 + c)e~3 b3_6c . Together with the boundary condition 
lim c _ KX> / / (c) = 0,weget (jQ). 

(C) Let f(c) := l.h.s.f lA.4]) . Differentiating twice and applying (1 A. 1 [) 



we obtain f"(c) = Ai(6 2 + c)e~^ bS ~ bc . The boundary conditions are 
lim^oo /'(c) = = linic^oo f"(c). Integrating twice and shifting the inte- 
gration bounds to zero, we get (IA.4|) . 

(D) This can be derived by ( ]A.3j) by change of variable W := —Z. 

(E) This identity can be found by first using flA.8j) and then the two 



complex integrals are decoupled. For the integral over W one uses flA.lj) . 
while for the integral over Z one uses (IA.3I) . 



(F) One applies the identity (2.20) of [2T], which holds for b 2 < h, 



[ c/Ae- A(i ' 2 - fel) Ai(c/ 1 + A)Ai(c/ 2 + A) 



{d 2 -d l f . d x + d 2 . (61-62) 1 . . 
ex P ~77Z TT _ ( & i _ b V n 1 77; ( A - 9 ) 



Wb 1 - b 2 ) V 4(6i - 6 2 ) v ZJ 2 12 
and then we set d\ = b\ + c\, d 2 = b\ + c 2 . □ 

B Invertibility of 1 - P s K Ai P s 

Lemma B.l. For any fixed real numbers si, . . . , s m , 

(1 - PsKmPs)- 1 (B.l) 

exists. 

Proof of Lemma IB. 11 First of all, notice that 

det(l - P s KmPs) = p( f) {^(r fe ) - r 2 < s,} V (B.2) 

^ fe=i ' 
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Now, let S := min{si, . . . , s m }. Then, 

mM>v(f){A(r k )-T*<s}) 



'- 1 ' (B.3) 

> P ( su P (.A(t) - r 2 ) <s]= F GOE (S) > 0. 



T& 



The equality with the GOE Tracy- Widom distribution, -Fgoe, is proven in 
Corollary 1.3 of |2T], while the strict inequality follows from the monotonic- 
ity of the distribution function and the large —S asymptotics (see for exam- 
pie, [3]). 

Moreover, as shown in Section 2.2 of |21j . P s KaiP s is trace-class (the shift 
by r| is irrelevant for that property, since it holds for any Sfc's). The proof 
ends by applying a known result on Fredholm determinant, see e.g. Theorem 
XIII. 105 (b) in [30J: let A be a trace-class operator, then 

det(l + A) ^ 1 + A is invertible. (B.4) 

□ 



C Trace-class properties of the kernel 

In this section, we discuss the basic properties of the operators K and K in 
Section [2j 

Proposition C.l. 

(i) The operators P U KP U and P U KP U are bounded in L 2 ({1, . . . ,m} x R). 
(ii) Let a,b E (0, |). Fix constants ati, . . . , a m such that 

- - < Qfi < a 2 < ■ ■ ■ < a m < (C.l) 

T/ien P U K P u defined in (12.21jl a trace-class operator on 
L 2 ({1, . . . ,m} x R). On the other hand, P u K con] P u defined in fT2TT8l) zs 
a trace-class in the same space if the constants satisfy more restricted 
condition 

— a < a.\ < a 2 < ■ ■ ■ < a m < b. (C.2) 
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Proof of Proposition IC.ll We only prove (ii). The proof of (i) follows 
easily by suitably modifying the analysis for part (ii). 

We start with P u K c ° ni P u . Since the set of trace-class operators is a linear 
space, it is enough to prove that P Ui MiKijM~ l P Uj is a trace-class operator 

in L 2 (R) for each As Kij = K^j — Vij (see (12. 13ft ). we prove that each 
of the two operators is trace-class. 

(a) Operator : It is easy to check that the operator can be re-expressed 
as (see (OS]) ) 

P Ui (x)e- a ^k iJ (x,y)e a ^P U] (y) = [ dzL(x, z)R{z,y) (C.3) 

Jr 



where 



L{x,z) = P Ui (x)e- a * x -± & dwU™>~ w(x+z) Po(z), 



2m 



R(z,y) = e «i»P tt .( y )-L I ^r__p o( -). 

J 1 -1/2 



1/2 r + , (C4) 



By choosing the integration paths close enough to 1/2, resp. —1/2, a simple 
estimate shows that there exists < 5 < min{ai + 1/2, 1/2 — a m } so that 

\L(x,z)\ < c e -(V^> e -WM 2l[ ^ 0]) 



From this it follows immediately that L and R are Hilbert-Schmidt operators 
on L 2 (M). So is the conjugated kernel of K it j trace-class. 

(b) Operator V$j: As Vij = for i < j, assume that i > j. 
In this case, we have 

At ^U-tj >1. (C.6) 

From (12.141) . by plugging in fa and <f>j and by changing the contour from ilR 
to R, we obtain 

1 r°° p -iz(.x-y) i roc 

- s L (in^iF^P ^ = s L e mdz (C7) 



where 



(i + lZj (f ~~ 12 
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Hence 

V iJ (x,y) = f(x-y), (C.9) 

the Fourier transform of /. 

Note that f(z) = g(z)h(z) where 

»W = (TTuP' h{z) = (F^' (C10) 

As g G L 2 (R) so then is g G L 2 (IR). Moreover, as 

1 f e~ izs 



dz (C.ll) 



27ry H (|+i^* 

is well-defined pointwise except at s = (since it is conditionally convergent), 
g(s) equals (IC.llj) almost everywhere. Furthermore, we can compute (1C.11I) 
explicitly and obtain 

g(s) = {At _ lV eS/2l i*<°]- ( c - 12 ) 

Similarly, 

Ks) = ^ e-'/ 2 l [8>01 . (C.13) 

We have 

(g*hy(s)=g(s)h(g) = f(s). (C.14) 

Hence 

Vi,j(x, y) = f(x -y) = (g *h)(x- y) 

•{x — y — s)h(s)ds = / g{x — s)h(s — y)ds 



(C.15) 



Here all the steps make sense since g,h G L 2 (R) fl L x { 

Note aj > ctj for i > j, and 1 ± (ctj + otj) > 0. By using flC15j) . we see 
that 

P^MiVijM-'Pu. = {P^M&N-^iNUM^P^) (C.16) 
where Q and H are the operators with kernel 

} , , , (C17) 

= ( t . _ t . _ X ), e 2 1 [^-2/>o]> 
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and TV is the multiplication operator 

N(x) = e^( '+ a >. 

Then 

[I dxdyKP^M&N-^&y^dydx 
J Jm? 



(C.18) 



dx dy 



y>x>Ui 



ai+a 3 )y 



(M-tj-i)\ 

Changing the variable y by s := y — x, the above equals 
1 



(C.19) 



dxe -(^- a] )x / dss 2{u-t^i) e ~(i^ ai ~ Qj )y (Q20) 



which is finite. Hence P Ui MiQN 1 is a Hilbert-Schmidt operator. Similarly, 
// dxdy\{NUMj l P U3 ){x,y)\ 2 

J JR 2 



dx dy 



x>y>Uj 
1 



e 2 



(x-y) 



{u-tj-iy. 

POO 

dxe- {a *- a > )x / rf ss 2 ^-*i- 1 ) e -( 1 +^+ a i)^ 



(C.21) 



((t,,-t,-i)!)2; u 

is finite, and hence Nl-LMJ 1 P u . is a Hilbert-Schmidt operator. Therefore, 
MiPuVijPujMj 1 is a trace-class operator. 

Now we consider i^ con j. From f !2. 16j) . 

P Ui K™pP U] = P Ul KpP U] + (a + h)P„ M,f, ® g j M- 1 P u 

But by changing the contour in (I2.12p . 



(C.22) 



fi{x) = <fri(a)e ax + — & dw 



2vri 



we 



r i/2 



a — w 



(C.23) 



Hence there is a constant C > such that < Ce~ ax for some constant 

C for x > Therefore, Mi(x)fi(x) G L 2 ((wj,oo)) if a, > —a. Similarly, 
gj(y)Mj(y) G L 2 ((w,-, oo)) if aij < 6. Being a product of two Hilbert-Schmidt 
operator, P u% Mifi <g> gjM~ 1 P U] is trace-class if (lC.2j) holds. This completes 
the proof. □ 
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D Gaussian fluctuations 

Consider the directed percolation model defined in fll.ip . but look away from 
the characteristic line. We set 

x=-l-N, y = -^—N. (D.l) 
1+7 1+7 

Denote Q(a, b) = G((a, b), (x, y)) the passage time from (a, b) to (x, y). Then 

G(x, y) = max{Q(0, 1), Q(l, 0)}. (D.2) 

In Section 6 of j2] (see also [27]) a last passage percolation with one source 
only is considered. Nevertheless, the arguments given therein can readily be 
used to prove the following: if 7 > 7 C = p 2 / (1 — p) 2 

lim P(Q(1, 0) < ciJV + c 2 sN 1/2 ) = -^L [ dx e~ x2/2 = $(s), (D.3) 

N ~>°° V27T J-oc 

where c x = + ^y) and c 2 = J^L-^- IX^. If 7 < 7c then 

there exists a constant c' 2 (see such that 

lim P(Q(1, 0) < c^iV + c'^iV 1 / 3 ) = F GUE ( S ), (D.4) 

iV— >oo 

where c' x = y^;(l + -^) 2 an d -^gue is the GUE Tracy- Widom distribution. By 

symmetry, if 7 < 7c , setting 61 = i^(iz^+^) and 6 2 = ^^^-^2 - ^2, 
it holds 

lim P(Q(0, 1) < &iiV + 6 2 siV 1/2 ) = $(s) (D.5) 

A 1 "— >oo 



and for 7 > 7 C 



lim P(Q(0, 1) < c\N + c' 2 siV 2 / 3 ) = F GUE (s). (D.6) 



Consider now the case 7 > 7 C . Then since c[ < C\ we have 

lim P(Q(0, 1) < Cl N + casA^ 1 / 2 ) = 1. (D.7) 

N— >oo 

With the notation d := ciiV + c 2 siV 1//2 and (1D.2j) we obtain 

P(Q(1,0) < d)>¥{G{x,y) < d) 

= P(Q(0,1) < dnQ(l,0) < d) 
= P(Q(0, 1) < d) - P(Q(0, 1) < d n Q(l, 0) > d) 
> P(Q(0, 1) < d) - 1 + P(Q(1, 0) < d). 
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(D.8) 



We take the limit iV — > oo on both sides and by f]D.3|) and (ID .70 we get 



lim F(G(x, y) < c x N + c 2 sN 1/2 ) = (D.9) 

N— >oo 

Similarly, for the case 7 < 7 C , one shows in the same way that 

lim F(G(x, y) < hN + b 2 sN 1/2 ) = $(s). (D.10) 

N— >oo 
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